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Abstract 

An increasing number of articles during the last five 
years have suggested that the atmosphere of the North 
American Arctic is subject to high aerosol loads in winter 
and low aerosol loads in summer. The winter maximum is said 
to be due to pollution probably originating in Eurasia. 

To confirm or reject this hypothesis the aerosol 
optical depth was calculated at five, and the aerosol 
single- scattering albedo at one station in the Canadian 
Arctic. Use was made of hourly radiation data and synoptic 
observations. 

The aerosol optical depth was calculated from a direct 
radiation model at Resolute and an empirical clear sky 
global radiation model suggested by Davies and Hay (1980) at 
the other four stations. 

The single-scattering albedo was calculated at Resolute 
using the above-mentioned clear-sky global radiation model, 
a two-stream approximation in a one-layer atmosphere, and a 
two-stream approximation in a_ two-layer atmosphere with 
variable aerosol and water vapor distributions. 

It turned out that the total amount and the seasonal 
march of the aerosol optical depth is similar at all _ five 
stations. The most outstanding feature is a short period of 
high aerosol turbidity in spring, followed by an abrupt 
decrease to low summer values. The southern and eastern 
stations have a second maximum in late summer. The few 


winter values suggest that winter is a period of rather low 
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aerosol turbidity. 

The aerosol single-scattering albedo was found to be 
about 0.8. As this result has a large error of possibly as 
much as +0.2 it cannot contribute to answering the question 
of the origin of the aerosol. The limitation of these models 
to calculate the aerosol single-scattering albedo is mainly 


due to a 5%-error in the radiation measurements. 
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introduc Fron 

One of the threats facing mankind during the late twentieth 
century is inadvertent climatic change. It is believed to 
have its origin mainly in two processes: the increase of 
both carbon dioxide and aerosol concentration on a 
wor ld-wide scale (Rasool and Schneider, 1971). Carbon 
dioxide is of crucial importance to the’ transfer of 
terrestrial radiation, while aerosol influences the solar 
spectrum. 

Carbon dioxide concentration is Known to be largely 
uniform in the atmosphere. This is not true for atmospheric 
aerosols which, due to their comparatively short residence 
times in the lower atmosphere, are concentrated in_ their 
natural and anthropogenic source areas. 

Thus, until recently, remote areas were said to be 
subject to a background aerosol load (Porch and Radke, 1970) 
which reflects the equilibrium state between natural aerosol 
production and destruction. It was not before the early 
1970s when measurements indicated that, even for a_ site as 
remote as northern Alaska, the atmospheric turbidity due to 
aerosol is, at least seasonally, comparable to non-urban 
points in midlatitudes (Shaw and Wendler, 1972; Radke et 
al., 1976). During the 1980s a multitude of research papers 
are revealing the spatial and seasonal patterns in the 
western Arctic and the chemical composition of the 
atmospheric particulates. The next chapter details’ the 


assumption that, in the western Arctic, the atmospheric 
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turbidity due to aerosol, in short aerosol turbidity, 
undergoes an annual cycle with a marked maximum in winter 
and spring, and originates from the midlatitudes. 

At the present stage of the discussion, the results are 
preliminary rather than definite. Porch and MacCracken 
(1982), for example, complain about ‘ the significant 
uncertainty in the optical constants of the arctic aerosol’. 
Therefore, more research is required into the areal extent 
as well as into the composition of the arctic aerosol. This 
is closely related to the question of the origin of the 
aerosol. 

From the point of view of climatology, a period of at 
least several years has’ to be investigated to reach more 
general conclusions. This is only possible if use is made of 
data obtained on a routine basis. These data are usually not 
tailored for the use in studies focussed on a special topic 
but rather for a general climatological overview. Therefore, 
the possibilities of retrieving sophisticated results, 
especially in the field of chemical composition of aerosols, 
are inevitably limited. 

The Atmospheric Environment Service, Canada, publishes 
hourly data sufficient to calculate the aerosol turbidity at 
about one dozen of arctic stations. For only one such 
station, Resolute (74.7°N, 95.0°W, see Figure 1) can the 
published data be used to calculate the single-scattering 
albedo, the interpretation of which allows conclusions as to 


whether the aerosol is heavily influenced by man-made 
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Figure 1: Location of some place names mentioned in this thesis 
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Dot lution: 

Theapurpese ofaathis mthesisetis = to cadculatessfor, a 
three-year period (uJanuary 1978 to December 1980) on an 
hourly basis the aerosol optical depth for Resolute, Alert 
(82.5°N, 62.3°W), Inuvik (68.3°N, 133.5°W), Baker Lake 
(64.3°N, 96.0°W), and Frobisher Bay (63.7°N,68.6°W), and the 
single-scattering albedo for Resolute using various 


empirical and theoretical radiative transfer approaches. 
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2. Literature Review 
Solar radiation data collected by the Atmospheric 
Environment Service consist of hourly global radiation 
readings covering the entire spectrum and, in the case of 
Resolute, of additional diffuse and reflected radiation 
readings of the entire spectrum. With these data only’ the 
aerosol optical depth can be calculated, which indicates the 
degree of turbidity of an atmosphere due to _ suspended 
particulate matter. For Resolute, the single-scattering 
albedo can be inferred from various radiation models. It 
indicates the contribution of aerosol absorption to aerosol 
attenuation. As natural aerosols are usual ly almost 
non-absorbing and man-made particulates may be highly 
absorbing, it can, as a first approach, be interpreted as a 


measure of the anthropogenic contribution to turbidity. 


2.1 Aerosol optical depth 
The attenuation of direct solar radiation can be 


described symbolically by 


I=T, 1,cos@ 2eule) 


where I, represents the unattenuated direct beam flux 
density (solar constant), I is the direct beam flux density 
on a horizontal surface at ground, O is the solar zenith 


angle and T, is the total transmissivity. In the solar 
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spectrum, T, can be approximated by 
og eta ee wheal 


with 1, being the transmissivity after ozone absorption, T, 
the Peace eery ity after Rayleigh scattering, Ty the 
transmissivity after water vapor absorption and T, the 
transmissivity after aerosol extinction. 


The latter can be represented by 
leosileuiee parecy 
with Ta. being the transmissivity after aerosol absorption 


and T.as the transmissivity after aerosol scattering. 


The transmissivities are usually written as 


T. =exp(-Tm) (254) 
T,, =exp (-T,,m) Gee) 
Tas =@Xp (-T,,m) (2.6) 


where = tw, t.. and t,, are the aerosol optical depth, the 
aerosol optical depth due to absorption, and the aerosol 
optical depth due to scattering, respectively. m is the 


relative optical air mass which is approximately 
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mz(cos@) (24m) 


As can be seen from equation 2.4, T, is independent of 
the solar position and is therefore an absolute measure of 
aeroso] turbidity. Thus, it can be used for comparison 
between different locations and different seasons. 

(2.2, can only be measured if the solar disk is 
unobstructed by clouds. But in many instances, the aerosol 
turbidity at ground level has been measured with integrating 
nephelometers which use an artificial light source and a 
small sample air volume. (see e.g. Charlson et al., 1967). 
This device allows the direct measurement of the aerosol 


volume scattering coefficient (Jas which is defined as 
d1/1=-f3,,dx (2,8) 


where dI/I is the fraction of radiation scattered out of a 
beam by the aerosol on a light path of length dx. /3., is a 
function of the particle number-size distribution and 
therefore generally a function of space. Introducing zZ as 
the vertical coordinate in the atmosphere, the relation 


between T,, and fas is then 


Tas= [Bae z)dz (2.9) 
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Aircraft measurements indicate that for the lower 5000m of 


the atmosphere /3.s(Z) can be approximated by 
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fS,,(z}=f,.(0)exp(-z/H, ) PLete: 


with H, being 1000m over central Europe (Penndorf, 1954), 
1400m over Alaska (Shaw, 1975), and 2000m over Nebraska 
(Blifford and Ringer, 1969). For a model standard 
atmosphere, Elterman (1964) suggests 1200m. From equation 


2.9 and 2.10 follows that 
G, 2Hss,,( OW (2141) 
or as a very first approximation 


T= 1200/3; 0) (Daas 23) 


with /3,,(0) expressed in m7’. 
In Chapter 2.2 it will be shown that 7. is usually not 
negligible compared with Ta; . For very preliminary 


estimations it may be assumed that on the global average 
ee) 25 )ce, (12.513) 


so that from equations 2.12 and 2.13: 


Ce =, At = 500 /08(.0) (2.14) 
Table 1 presents an incomplete list of some measurements of 


T. for various regions of the world excluding the Arctic 
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made during the last few years. Values in brackets indicate 
calculations made from (Sas (0) -readings of integrating 
nephelometers using equation 2.14. 

At this juncture it has to be pointed out that t,,and 


fsa, are also functions of wavelength 


Tas (A) 10.08; aClnatirsin 
Toate) ees (As) (+) eae 
Angstrém (1961) found that the wavelength exponent, «, is 


1.3+0.2 but may be as low as 0.5 after volcanic eruptions 
and forest fires (see, e.g., Ryznar et al., 1981). Many 
subsequent studies confirmed Angstrom’ s findings. 

Unfortunately, t. and as 0) are some t imes being 
determined for A=0.5um and sometimes for the entire solar 
spectrum. For Table 1, no attempts have been made to 
calculate t, for the entire solar spectrum from the values 
given for A=0.5mm. But with regard to the uncertain 
assumptions made so far, measurements made at around 0.Sum 
may be taken to be similar to those for the entire’ solar 
spectrum. 

Table 2 presents various values of T, for Barrow, 
Alaska, and Table 3 for other stations in the western 
Arctic. Values in brackets are again estimated from 
integrating nephelometer measurements of the scattering 
coefficient at ground level, (as(-0) , using equation 2.14. It 


is obvious that the values in brackets are approximately by 


a factor of 3 lower than the values gained from direct 
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Barrow, Alaska (71.3°N, 156.6°W) 


a Aum 


(0.05-0. | (0.05-0.09 Radke et al., | Radke etal, 1976 


Shaw, 1975 
C105" 022 
(0.02) southerly winds 
Peterson et al, 1980 
(0.05) northerly winds 
197771978: 
Jan. (0.02) 
Feb. (0.04);(0.02) 
Mar. (0.02);(0.02) 
Apr. 7 (0.02);(0.02) 
May (0.0 1) 
Jun. (0.002);(0.003) 
Jul. (0.009) 
Aug. (0.006);(0.005) 
Sep. (0.004);(0.004) 
Oe (0.008);(0.006) 
Nov. (0.0 1);(0.0 1) 
Dec. (0. 1);(0.009) 


Bodhaine et al., 1981 


Us =0:2 

Op —Oi2 

UO5—0 41 Shaw, 1981 
G:02—0.05 

0.02-0.04 


Table 2: Aerosol optical depth T,. for Barrow, Alaska Values in brackets were calculated 
from aerosol scattering coefficients using equation 2. 14. 
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Brooks Range, Alaska (69.3°N, 143.8°W), elevation 1740m 


1970-/4: 

1st half May 

2nd half Jun. Shaw and Wendler, 1972 
2nd half Jul. 


Mould Bay, Canada (76.2°N, 119. 3°wW) 


1979; 1980: 

Jan. (0.05) 

Feb. (0.06) 

Apr. (0.03);(0.04) 

May (0.02);(0.03) 

Jun. (0.005) 

Jul. (0.003) Barrie et al, 1981 


Aug. (0.005) 
Sep. (0.0 1) 
Oct. (0.02) 
Nov. (0.02) 
Dec. (0.04) 


continued on next page 
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iosa 


Barrie et al., 


aga Svalbard (79°N, 12°E) 


Apr.,May Bit th ee (0.02) Perens 1980 


Table 3: Aerosol optical depth T. for various arctic stations, except Barrow, Alaska 
Values in brackets were calculated from aerosol scattering coefficients using equation 
2.14. The asterisk indicates values which had to be calculated from Angstrém turbidity 
coefficients (Angstrém, 1961) assuming a wavelength exponent of 1.3 
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observations of the direct component of the solar radiation 
flux density. But for other regions of the world (Table 1) 
the values in brackets seem to be compatible with the direct 
measurements of 7@,. 

Therefore, it must be concluded that equation 2.14 is 
Notre appd icabileseforenaretictistationss adihisreisesurprising 
because equation 2.10 was found to be correct for Barrow, 
Alaska, with a scale height H, of 1400m, as mentioned 
before. Also, the errors incurred by the use of equation 
2.13 can certainly not explain the difference. Chapter 4 
will give as possible explanation that the scale height Ha 
is, contrary to the findings of Shaw (1975), probably 
several times greater. 

Thus, integrating nephelometer measurements cannot be 
used to determine the aerosol optical depth in the Arctic. 
But only the aerosol optical depth can be used as a measure 
of the aerosol turbidity of the entire atmosphere, because 
determination of the aerosol optical depth is an automatic 
integration of the local turbidity parameter /3,,(z) from the 
bottom to the top of the atmosphere, as can be seen from 
equation 2.9. 

Comparison of Tables 2 and 3 with Table 1 reveals two 
striking features: Firstly, the average optical depth in the 
Arctic is distinctly below the average of the midlatitudes. 
This statement is based on the interpretation of the 
directly measured aerosol optical depths (values without 


brackets). Secondly, there is a marked seasonal march of the 
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aerosol optical depth, with the lowest values in summer, a 
continuous increase in fall and winter, and a maximum in 
spring. A dramatic decrease between May and June is the most 
characteristic feature of the seasonal pattern. 

The seasonal turbidity pattern is not identical every 
year, as can be seen from the extreme difference at Barrow 
between December 1977 and December 1978. Differences on a 
day-to-day basis seem to be attributable to different air 
masses, as can be seen from aerosol turbidity differences 
between days with southerly and northerly winds at Barrow in 
winter. Similar observations were made by Unsworth = and 
Monteith (1972) for the British Isles and by Fitzgerald 


(1980) on a voyage across the North Atlantic. 


2.2 Single-scattering albedo 
The single-scattering albedo ™, for an aerosol is 
defined by 


ipl | (2.16) 


Bas + Baa 
with the aerosol volume scattering coefficient, /3,;, as 
defined in equation 2.8, and the aerosol volume absorption 
coefficient, ,,, defined in an analogous manner. 
In many instances, however, not the single-scattering 
albedo ~ but the imaginary part n, of the complex index of 


refraction 
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n=n,-in,z (2.17) 


is being determined by laboratory experiments with collected 
aerosols. 


A functional relationship between n and w, can be 


2 
found for a given aerosol number-size distribution, a given 
real part mn, of the complex refractive! index vof the aerosol’, 
and a given wavelength A. | 

The derivation of this functional relationship 
necessitates a complete solution of the Mie scattering 
problem. Besides the approach by Hansen and Travis (1974), 
only Twomey (1977) has calculated this relationship for a 
Junge distribution of particles between 0.01 and 10um, with 
a real part of n,=1.5 at A=0.55um. The relationship is shown 
in Figure 2. Twitty and Weinman (1980) calculated w, at 
A=0.5um for carbonaceous material (n=1.8-0.5i) and two 
non- Junge size distributions which are supposed ly 
representative of urban aerosols. 

Simultaneous measurements of ™4 and n, were- only 
reported by Chin-I Lin et al. (1973) for the New York City 
aerosol and by Roessler and Faxvog (1979) for acetylene soot 
produced in the laboratory. These latter measurements 
indicate that Twomey’s calculations can be used to estimate 


WZ from on with an error of less than 0.1. In fact, Junge 


72. 
aerosol number-size distributions have been found to be good 
approximations for many locations (cf. recent findings for 


Denver (Willeke and Whitby, 1975) and Mildura, Australia 
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(Gras and Michael, 1979)). 

Table 4 presents values for the single-scattering 
albedo w, for various aerosol substances. Table 5 presents 
ZW -values for atmospheric aerosol samples taken at various 
locations. In Tables 4 and 5, values in brackets’ indicate 
calculations From Nn, “measurements using Twomey’s 
relationship (Figure 2). They have to be treated with care. 
Measurements were being made at different wavelengths. But, 
as can be seen from Twitty and Weinman (1971), the 
single-scattering albedos for two different non-Junge size 
distributions are only slightly dependent on wavelength for 
A < 3um. 

It is obvious from Table 4 that the dominant natural 
aerosols, viz. sulfates, sea salts, and crustal particles, 
all have single-scattering albedos close to unity. OQn_ the 
other hand, products obi combustion show very low 
single-scattering albedos around 0.10 to 0.15. In other 
words, products of combustion are highly absorbing, whereas 
natural aerosols are almost. pure scatterers. Sadler and 
Charlson (1981) found a high correlation between the aerosol 
volume absorption coefficient /3,, and the carbon content of 
an aerosol in Washington state. 

But, with the exception of forest fires, products of 
combustion are almost exclusively man-made. Thus, the 
single-scattering albedo should be a tool to determine to 


what degree an aerosol is influenced by man’s activities. 
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0999 


0.999 visible Toon and Pollack, 1976 


0.98-0.96 visible 


acetylene soot, 


diesel exhaust 


(0. 12) O51 
Roessler and Faxvog, 1980a 
(0.12) Oe 


Table 4: Single scattering albedos we for various aerosol substances. Values in brackets 
are calculated from measurements of the imaginary part of the refractive index using the 
functional relation of Figure 2. 
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20 
US: 


OS0-0'S5 
0: 73=0:87 
0.91-0.94 
0.68-0.74 
0.89 
(0.60) 


7 urban-industrial sites 


7 urban-residential sites 


Waggoner et al., 1981 


3 remote sites 


New York City Chin-| Lin et al, 1973 


St. Louis Method and Carlson, 1982 


fly ash from forest Grams et al., 1972 


fires in the U.S. 


(0.70) 
(0.90) 
(0.90) 
(0.90) 


Tucson, Arizona King, 1979 


Big Spring, Texas Grams et al., 1974 


DeLuisi et al. 1976 


Big Spring, Tex., and Blythe, Cal. 


Lindberg and Laude, 1974 
Lindberg, 1975 


desert of southern New Mexico 


(0.90) 
0.80 
Oke}, 
(0.90)—(0.95) 
(0.85) 


Gainesville, Florida Ward et al., 1973 


southern England Roach, 1961 


Shetland Islands, Great Britain Forbes and Hamilton, 1971 


Irish west coast 


Fischer, 1973 


Jungfraujoch (3560m), 
(Switzerland) 


Fischer, 1970 
Eiden, 1966 
Fischer, 1973 


(0.65)—(0.85) 
(0.45)—(0.85) 
(0.65) 


industr. aerosols, Germany 


Mainz, Germany 


Mainz, Germany 


North Atlantic: 


Saharan dust Carlson and Caverly, 1977 


Saharan dust 


Patterson et al., 1977 


sea salt aerosol 


(0.65) Fischer, 1973 
(0.40) 
(0.70) 


(0.95)-(1.00) 


‘push’, near Tsumeb, Namibia 


Tel Aviv, Israel 


Levin et al., 1979 


Negev desert, Israel 


desert of Iran and Pakistan 0.5-1.1| Otterman et al., 1982 


Table 5: Single scattering albedos Wo for various locations . Values in brackets are 
calculated from measurements of the imaginary part of the refractive index using the 
functional relationship of Figure 2 and rounded off to the nearest multiple of 0.05. esis. 
means the entire solar spectrum. 
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Table 5 shows that this prediction holds approximately. 
The lowest values of ©, that have been measured are around 
0.40 to 0.65 and can be found near to urban-industrial 
sites. At urban-residential sites and cities without major 
pollution sources W, may be expected to range between 0.75 
and 0.90. Values around 0.85 to 0.90 are already typical for 
rather remote sites. But single-scattering albedos above 
0.90 can mainly be found in pollution-free and mainly 
maritime air masses, e.g., on the western coast of Ireland. 

Some values from Table 5 need a comment: The _ forest 
fire fly ashes’ show values as low as 0.60, hence they are 
indistinguishable from man-made combustion. products at 
heavily polluted sites. The value of 0.59 for the Shetland 
Islands, if reliable, indicates’ that single-scattering 
albedos typical of heavy pollution can, in fact, be found at 
a distance as far as 1000kKm from the major. pollution 
sources. The low value of 0.65 in the semi-desert near 
Tsumeb, Namibia, was measured during the dry season and is 
due to organic matter, viz. dry grass, being ground by wind 
erosion to aerosol-sized remnants and lifted into the air 
(Fischer, 1973). This organic material, although not burnt, 
behaves optically like ashes. This’ may also be an 
explanation for the low single-scattering albedos measured 
at other desert locations. 

It should be pointed out that neither direct nor 
indirect measurements of the single-scattering albedo are 


available for any arctic station. 
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The single-scattering albedo is a quantity of crucial 
importance in global climatic studies. Until about twenty 
years ago, it has been tacitly assumed that aerosols are 
virtually non-absorbing, and hence the _ single-scattering 
albedos are close to unity. 

It was not before the sixties when Roach (1961), 
Robinson (1962,1966), McCormick and Ludwig (1967), and 
Charlson and Pilat (1969) drew attention to the absorbing 
nature of aerosols. In the case of aerosols which both 
scatter and absorb, any increase of the total load also 
results in an increase of absorption in the atmosphere. 

Several subsequent studies investigated whether a 
higher aerosol load leads to falling or rising temperatures 
of the surface-atmosphere system in the case of a_ cloudless 
atmosphere without heat advection: 

From the model given by Ensor et al. (1971) a critical 
single-scattering albedo ©%.- can be calculated which is a 
function of the surface albedo, a, and of the ratio of 
aerosol backscattering to total scattering which is a 
function of the solar zenith angle. This ratio will be dealt 
with in detail in Chapter 4. 

If the single-scattering albedo of the aerosol w, 
exceeds 4 .-, an increasing aerosol concentration entails 
falling temperatures of the surface-atmosphere system and 
vice versa. For a solar zenith angle of approximately 1H0Pe 
and surface albedos of 0.2 and 0.8, respectively, the 
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Wo cr = (e248) 


Atwater (1970), in a similar approach, comes up with a 
result from which follows that ™.- is also a function of the 
tat1o;, R, of global radiation at the ground to global 
radiation at the top of the atmosphere. Therefore, for a 
given solar zenith angle (e.g. 70°), and given combinations 
of aandR (e.g., at Resolute, typical values for cloudless 
hours are: a=0.2 together with R=0.7, and a=0.8 together 


with R=0.8, each for zenith angles of @=70° ) “ocr can be 


calculated to be 


We,er = (e299 


As, by definition, 0s%<1, the very last result means’ that 
any aerosol increase leads to a warming of the 
sur face-atmosphere system. It turns out that the value _ for 
Wye from equation 2.19 becomes identical with the value for 
Weer from equation 2.18 for the rather unrealistic assumption 
of R>0, i.e., an atmosphere of infinite optical thickness. 
Coakley and Chylek (1975) present a result which can be 
found to be also aé_e function of R and which is exact ly 


correct only for R71, .i.e., a completely transparent 
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atmosphere. 

From Twomey (1977,p.285f.) a critical single-scattering 
albedo can be found which is independent of RikFor’ ale solar 
zenith angle of 70° and albedos of 0.25 and 0.75, 


respectively, the following holds: 


0.80 for 1a=0725 
Wo cr = (12322'0;) 
0.99 forwacs0 2/5 


Yamamoto and Tanaka (1972) find that an imaginary part 
of the refractive index of 0.05 (corresponding to #,=0.6, 
according to Figure 2) is approximately the critical value 
for low albedos (0.05-0.15) and various zenith angles. 

In a very extensive study containing a comp lete 
solution to the radiation transfer problem, Herman and 
Browning (1975) presented a result indicating that, when 
a=0.8, warming occurs in any case, whereas for a=0.1 and a 
solar zenith angle of 55° the critical imaginary part is 
0.018 (corresponding to 4,=0.79). 

There are several more attempts to solve the problem of 
climatic change with increasing aerosol load. But’ the 
results of these studies are not reliable because important 
terms were forgotten right at the outset of the 
calculations: Mitchell (1971), Rasool and Schneider (1971), 
Chylek and Coakley (1974), and Russell and Grams (1975) all 


forgot to take Rayleigh scattering into account. This is 
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under no circumstances justifiable, as Rayleigh scattering 
is usually at least of the same order of magnitude as 
aerosol attenuation. Neumann and Cohen (1972) do not 
consider’ reflection» of “global radiation fromthe ground; 
thus their model is not applicable for snow covered 
surfaces. Finally, Barrett (1971) assumes that aerosols are 
not absorbing, which is an invalid approach. 

So, it can be concluded that, for the high albedos 
typical for arctic snowfields, the addition of almost any 
aerosol leads to a warming of the surface-atmosphere system, 
whereas for the low albedos of the tundra in summer this 
question cannot be answered unambiguously. 

A more detailed study of the influence of aerosol on 
the terrestrial and atmospheric temperatures, specifically 
for an arctic environment, was made by Shaw and Stamnes 
(1980). It concludes that, after a layer of aerosol with 
typical values for ©, and ~, is introduced into the model 
atmosphere, the shortwave radiation balance is altered such 
that the atmosphere is being heated and the surface being 


cooled. The longwave radiation balance is not affected. Ina 


very similar study, Grondin (1980) confirms these results. 


2.3 Origin of the arctic aerosol 
The extent and seasonal variation of the arctic aerosol 
turbidity was depicted in Section 2.1. This section tries to 


elucidate the question of its origin. This question has a 
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political dimension because statements about the sources are 
automatically statements about the ‘responsibility’ for the 
so-called ‘arctic haze’. 

At first, arctic aerosol turbidity was thought to be an 
autochthonous phenomenon due to _ ice crystal formation at 
very low temperatures (Shaw, 1975). But, shortly thereafter, 
Rahn et al. (1977) discovered from chemical analysis that 
true aerosol layers are responsible for occasionally high 
turbidities. 

A relatively simple but still powerful tool to 
determine whether an aerosol is mainly composed of the 
materials of the earth’s crust, or is composed of sea salt 
or is of a different composition, is the enrichment factor. 
For crustal material, the enrichment factor E, with respect 


to any element X is defined by 


2 Cx Al accosal 


bs - (X7AI) crust 


(29218) 
where X means’ the mass of an element X and Al the mass of 
aluminum for the aerosol sample and for the average of the 
earth’s crust, respectively. The closer E, approaches unity 
the more reliable is the statement that an aerosol sample is 
made up of crustal material. For sea salt the enrichment 
factor is defined similarly. 

Using equation 2.21 Rahn et al. (1977) found that the 
events of high aerosol turbidity at Barrow in spring are 


probably due to a crustal aerosol. 
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In a later study Rahn (1981b) defined the noncrustal 


concentration of an element X in an aerosol by 


Dicer Maner Ott Goth a ( X/A1 res A ere (2 P22 ) 


With vanadium and manganese as elements X, appreciable 
noncrustal sconcentrations swere\ found) in othe Arctic for 
winter (Rahn and McCaffrey, 1980). Noncrustal vanadium stems 
from combustion Of residual oil, whereas noncrustal 
manganese has various sources. 


Now a ratio Y can be defined by 


Y=Ming, er / Vawter (228) 


Y was found to be 2.0 for Eurasian aerosols, but 0.4 for 
North American aerosols. For the Arctic, Y-values greater 
than 1.0, sometimes even as high as 2.0, were reported. With 
regard to the fact that Y decreases with increasing age of 
an aerosol, Rahn (1981b) concludes that the arctic aerosol 
originates in Eurasia rather than in North America. 

Rahn and Heidam (1981) in a general review on arctic 
air chemistry between 1977 and 1980 report that graphitic 


/ 


carbon was. found ‘in abundance’. Rosen et al. (1981) prove 
the existence of graphite by Raman spectroscopy and by 
determining the wavelength dependence of absorptivity, as 
well as the solubility and oxidation temperature of 


collected aerosol samples. This finding, together with the 
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finding that the aerosol size spectrum is typical of highly 
aged aerosols, iS one more indication that the arctic 
aerosol is a pollution-derived aerosol from aé_ée distant 
source, at least during the periods of high turbidity. 

From measurements made at Barrow in 1979 Daisey et al. 
(1981) report that particulate organic matter is in March, 
four times as abundant as in August; the aerosol is reported 
as being dominated by fly ash in March and pollen as well as 
sea salt in August. But Barrie et al. (1981) in a study made 
on aerosol samples taken at Mould Bay (76.2°N, 119.3°W) and 
Igloolik (69.4°N, 81.8°W) conclude that the pollen originate 
exclusively in North American deciduous forests. This is not 
only true for September and October, but also for April and 
May which is the season with the heaviest turbidity. Thus, 
North America has to be considered as a possible source 
region, too. 

A different approach to this problem is the synoptical 
one (Rahn, 1981a): The flow pattern in the Arctic in winter 
is generally determined by the Siberian high, the Icelandic 
low, and the Aleutian low. Therefore, there are two major 
pathways for midlatitudinal air to the North American 
Anctictartromancuropem) cuttings acrosse themwArctic Ocean 
somewhere between the North Pole and the northern coast of 
Asia, reaching Alaska and northern Canada as a northerly or 
northwesterly wind, and from eastern Asia across the Pacific 
Ocean and along the Canadian west coast, reaching Alaska as 


a southerly wind. Air from eastern North America could reach 
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the Arctic only via Europe. 

But these pathways are not equivalent in their efficacy 
for bringing high pollution loads to the Arctic. This is so 
because precipitation is the controlling factor for aerosol 
residence times, due to the processes of rain-out and 
wash-out. Both trajectories, from eastern Asia to Alaska, 
and from eastern North America to Europe, go over stormy and 
rainy seas. But the trajectory from Europe via the Arctic 
Ocean is a cold one with no major precipitation (Rahn, 
1981a). Remembering that northerly winds bring greater 
turbidity to Barrow in winter’ than do_ southerly winds 
(Peterson et al., 1980) the synoptic situation supports the 
hypothesis of Europe as the pollution source. 

In summer, these flow patterns do not exist, and_ the 
air in the Arctic is ‘decoupled’ (Rahn and McCaffrey, 1980) 
from midlatitudinal air masses through the polar’ front. 
Simultaneously, precipitation reaches its maximum, and so do 
the scavenging processes. This may explain the marked 
seasonal march of arctic aerosol turbidity. 

Another approach to the problem of the origin of high 
aerosol turbidity is the reconstruction of actual 
trajectories on a given pressure surface. Rahn et al. (1977) 
by calculating back-trajectories on the 70kKPa surface from 
Barrow, Alaska, for April and May 1976, arrive at. the 
conclusion that, at least in spring, direct connections 
exist between high turbidity in northern Alaska and_ flows 


from the central Asian deserts of Takla Makan and Gobi to 
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Alaska. 

Heavy soil dust loads with origin in eastern central 
Asia have also been found at Enewetak Atoll] (11°N, 162°E) in 
April 1979 by Duce et al. (1980) and on Hawaii (19°N, 156°W) 
from late April to early May 1979 by Shaw (1980b) and Darzi 
and Winchester (1982). Duce et al. assert that ‘in China 
dust storm activity is greatest in the spring because of the 
combined effects of low rainfall, the increased occurence of 
high surface winds which are associated with cold fronts, 
and soil freshly plowed for planting’. 

In this context, it should be remembered that the 
maxima of turbidity in the Arctic occur in spring at about 
the same time when the central Asian dust storm activity 
occurs. The transarctic flow pattern from Europe to Alaska 
mentioned above can certainly not explain such a maximum. 

More trajectory analyses were performed on the 85kPa 
level for Mould Bay from December 1979 to January 1980 
(Barrie et al., 1981) and for Barrow from February 1975 to 
January 1980 (Miller, 1981). At Mould Bay, 50% of the 10-day 
back trajectories originated in polar regions, the remaining 
50% came from Siberia, Europe, and North America. At Barrow, 
60% of the 5-day back trajectories had their origins inside 
a circle of 2000km radius around Barrow, with a lower 
percentage in winter and a higher percentage in summer. For 
the trajectories from outside the 2000kKm radius there were 
two distinct source areas: the northern Pacific (southerly 


winds at Barrow), and the Kara- and Barents-Seas (northerly 
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winds at Barrow). The north Pacific flow has a maximum in 
January and a minimum in June, with a secondary August 
maximum. The Eurasian flow is steady between July and March, 
and has a distinct minimum between April and June. 

At this stage, no conclusion can be drawn as to the 
origin of high turbidities in the North American Arctic. 
There are three possible source areas: Europe, North 
America, and China. Evidence is availabie in favour of each 
of them. Unfortunately, almost all investigations have been 
concentrated on Barrow, Alaska. Investigations from Svalbard 
(79° N,12° E)(Heintzenberg, 1980; Heintzenberg et al., 1981) 
are so far not comprehensive. With the exception of a_ study 
by Polavarapu (1978) which yields doubtful results, only one 
investigation on arctic turbidity has so far been done _ in 
the vast territories of the Canadian Arctic (Barrie et al., 
1981). 

The purpose of this thesis is to fill some of this gap 
with investigations on the aerosol optical depth at five 
locations and the single-scattering albedo at one location 
in the Canadian Arctic. This is a necessary contribution to 
the Knowledge of the spatial and seasonal patterns of the 
occurence of aerosol turbidity. Thus, it may be helpful for 


answering the question of its origin. 
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3. Aerosol optical depth at Resolute 


3.1 Direct solar radiation model 
In Chapter 2.1 the simplest model for the transfer of 
direct solar radiation was introduced (see equations 2.1 and 


cae 

Bales: (atecos © oes) 
or abbreviated 

I=JTT, 1, cos hers 


This multiplicative approach is, strictly speaking, only 
applicable if all transmissivities are independent of 
wavelength. This is certainly not true for the transfer of 
solar radiation in the earth’s' atmosphere, as ozone 
absorption can almost only be found in the ultraviolet 
spectral region, whereas water vapor absorption is 
restricted to the near infrared region. Rayleigh scattering 
obeys the wel1-Known X’-law, so its main contribution stems 
from short wavelengths. Even aerosol extinction is not 
wavelength independent: aerosol scattering was reported to 
obey approximately a Dean. The wavelength dependence of 


aerosol absorption is dependent on the aerosol material. 
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If the solar spectrum is divided into two equal parts 
at A=0. 73pm FCoUlSCOnme Ov 5 sb ecl > ee bycanabes assumed. shat 
ozone absorption and Rayleigh scattering is confined to the 
lower part of the spectrum, and water vapor absorption 
occurs only in the upper part thereof. Aerosol extinction is 
assumed to be equally distributed over the entire spectrum. 

Then, according to Paltridge and Platt (1976, p.121), 


equation 3.1 has to be written as 
Lette eas llecos:€ (358) 


where a, is the water vapor absorptivity, which is_ related 


to water vapor transmissivity T, by 
ageanehe Com4o 


This can be explained with an example: 

images al 0ncom On. (oly Ter 0 aS) soln =04 7 Se Oimame= UNO), 
and 1T, = 1.00. As extinction due to ozone absorption and 
Rayleigh scattering is supposed to be confined to the lower 
50% of the spectrum and water vapor absorption to the upper 
50% thereof, 30% (or: 2*0.15) of the lower half and 50% (or: 
2*0.25) of the upper half of the spectrum is subject to 
extinction,  “tre., 40% of the entire spectrum is 
extinguished, in other words, the total transmissivity is 
60%. Using the multiplicative approach (equation 3.1) would 


yield a total transmissivity of 0.85*0.75*1.00 = 0.64. But 
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equaltionigaseyields the iconrect solution: (0.85-0.25)+1.00 = 
0460. 
Remembering that the aerosol transmissivity T, can be 


expressed as 
T, =exp(-T,m) aoe 5) 


equation 3.3 can be solved for the aerosol optical depth 


epee sie 2 SP 
(To Te - au) I, cos 9 


T.=-= In| (3.6) 
This is the basic formula for calculating the aerosol 
optical depth from observations of direct solar radiation. 
This direct radiation model and the evaluation of its terms 


as described in the next section can be found in Davies and 


Hay (1980). 


3.2 Evaluation of the model 

To evaluate equation 3.6 seven quantities have to be 
Known : 
-the solar radiation flux density outside the atmosphere I, ; 
-the direct horizontal solar radiation flux density on the 
ground |; 
-the solar zenith angle @; 
-the optical air mass m; 


-the ozone transmissivity T, ; 
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“the Rdyleigwvuransmiussivitlyen., 


-the water vapor absorptivity a,. 


3.2.1 Solar radiation flux density outside the atmosphere 

The solar constant or solar radiation flux density 
outside the earth’s atmosphere at mean sun-earth distance is 
widely believed to be 1353+14Wm~* (Thekaekara and Drummond, 
171. This value is referenced to the International 
Pyrheliometric Scale (IPS 1956) (see e.g. Latimer, 1973, and 
Thekaekara, 1976). 

Recently, Jacobowitz et al. (1979) suggested 1368Wm 
based on Nimbus 6 satellite measurements. However, this 
value is referenced to a new cavity radiometer scale which 
is approximately 2% higher than the older IPS 1956 and has 
not yet found international recognition. As the measurements 
used for this thesis are referenced to IPS 1956 (Canada, 
Environment Canada etc., 1978ff.) the new value would have 
to be lowered by 2% to 1341Wm and would therefore be inside 
the error limits given by Thekaekara and Drummond (1971). 
Due to the uncertainties accompanied with the new scale, the 
old value of 1353Wm ~is being used throughout this thesis. 

The variable relative sun-earth distance, r, must be 
accounted for by multiplying the solar constant with a 
factor 1/r* to obtain the solar radiation flux density 
outside the atmosphere I,. Spencer (1971) calculated this 


factor as being 
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1/r>=1,000110+0.034221cos7*+0.001280sinr* 
0.000719c0s2V7+0.000077sin2y Len 74 


with 


V227d/365 (3.8) 


where d is the number of the day starting with d=0 on 


January ist and going to d=364 on December 31st. 


3.2.2 Direct solar radiation flux density on the ground 
Direct horizontal solar radiation is not being measured 
at any site in the Canadian Arctic. But at one arctic 
location, mamely Resolute, two pyranometers are working 
simultaneously: One measures global radiation from. the 
entire celestial dome, and the other measures diffuse sky 
radiation. This is achieved by shadow bands mounted in- such 
a way that the solar disk is obstructed. By subtraction of 
the two values direct solar radiation can be calculated. All 
pyranometers used in the Canadian Arctic are Eppley 
precision spectral pyranometers (PSP) (Latimer, 1980). 
Canadian measurements are published as hourly radiation 
flux density values starting at each full hour _ local 
apparent time (LAT) given in the units Mud/(m*hour). 
According to Latimer (1980), the root mean square (RMS) 
error of radiation readings is of the order of +5%. Under 


the best conditions it may be as low as +2%. 
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The implied error in ©, caused by the 1% inaccuracy of 
I, and the 5% inaccuracy of I will be discussed in detail 


later. 


3.2.3 Solar zenith angle 
The solar zenith angle O can be calculated from the 
well-Known relation (see, e.g., Kondratyev, 1969, p.342, or 


Liou, 980% p246) ° 
cos@=sin(v)sin(d)+cos(¥)cos(d )cos(h} G3). 9°) 


where “/ i's the latitude, b the solar declination, and h the 
hour angle. 
The declination can be expressed according to Spencer 


(1971) by 


5 =0.006918-0.399912cos7+0.070257sinv=0.006758cos20% 
0.000907sin27-0.002697cos37+0.001480sin37 (321.08) 


with ” as defined in equation 3.8. 
The hour angle h in degrees can easily be calculated 
from 


h=15* |12-LAT| is. 71 


where LAT means the local apparent time in hours. 
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In reality, the sun appears to be closer to the zenith 
than predicted from equation 3.9. This is due to the 
refraction in the atmosphere. If 9 denotes the astronomical 

/ 
zenith angle (equation 3.9) and GO denotes the zenith angle 
/ 
observed with refraction, the difference 9-9 can be written 
as (Cameron et al., 1963): 
nS gree aS) 
(ete Ts Veh ili sree ea (see) 
: “[=0 

WUEieE Dei nguthne abso lutemstlation sstemperature vidivided =by 


273K, r, being the station pressure divided by 101.3kPa, and 


? 
B.(@) being coefficients tabulated by Garfinkel (1944). 

The first and second columns of Table 6 present the 
ratio of the cosines of the real (O) and astronomical (0) 
zenith angles. For zenith angles of approximately 65° the 
difference is 1% and for zenith angles of approximately 85° 
the difference is as high as 10%. The error in ©, incurred 
by neglecting refraction can be calculated from equation 


Syiehe 


Te ONG. = | a 
AT,= s—*, 4 cos8= ——7(cosé cos0) (Sore 


The third and forth columns of Table 6 present the 
absolute changes in t. incurred by including refraction. The 
changes increase with increasing zenith angle and are about 
OT, =+0.01 at 85°. Because great zenith angles are quite 


common in the Arctic and the absolute values of Cn are, at 
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Table 6: Z=cos@'/cos9, with @' being the zenith angle with and 9 being the zenith angle 
without refraction, for a summer and a winter atmosphere (T=273K, and 243K on the 
ground, respectively). 

AT= T.'- Ta, with T.' being the aerosol optical depth with and T. being the aerosol 
optical depth without refraction, for a summer and a winter atmosphere. 

All caiculations were made for a station pressure of 101.3kPa 
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least seasonally, of the order of 0.05 and below, the 
refraction correction is not negligible. Refraction has so 
far been ignored in all radiation studies Known to the 
author of this thesis. 

The use of hourly readings of I in equation 3.6 


necessitates an integration of cos6(t) over one hour 


cos@= [cos@( t)dt/( thour) (Biat4 |} 

| hour 
This integration was avoided by using mid-hour values of 
cos 9(t), instead. The error incurred with this procedure is 
particularly low in the Arctic because cos Citi eedeesmynot 


change rapidly as the sun moves along the horizon. 


3.2.4 Optical air mass 

The optical air mass is a dimensionless quantity which 
is the ratio of the air-density weighted path length of a 
solar ray through the atmosphere to the air-density weighted 
path length of a ray from the zenith. 

Neglecting the curvature of _ the earth and the 
refraction of the atmosphere, the air mass m is simply given 


by 
meicos@)" bSest5:) 


This is a good approximation for small zenith angles 0. 
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For solar positions close to the horizon the effects of 
curvature and refraction are not negligible and have to be 
taken into account. Kasten (1966) calculated m as being 


- 1.253 


m=1/(cos0+0.1500(93.885-6) (3.16) 


where the zenith angle 0 has to be inserted in degrees. For 
6 =85° , m according to equation 3.16 is by 10% less than m 
according to equation 3.15, for @=80 the difference is 
still 3%, and for 9=70° the difference is already less than 
1%. 

To account for variable station pressures m may be 
multiplied by a pressure correction term yielding the 


corrected air mass m’ 


m’ =m(p/ 101.3) (321%) 


where the station pressure, p, has to be inserted in 
Kilopascals. This correction is important if C,-results from 
stations at different elevations are to be compared. 

From the definition of the optical air mass as an air- 
density weighted relative path length, it is obvious that 
its applicability is, strictly speaking, limited to Rayleigh 
scattering. This is so because the mixing ratios of ozone, 
water vapor, and aerosol are not constant with height. 
Seperate optical air masses should be calculated for every 


single hour, taking into account the actual vertical mixing 
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ratio profiles. Due to the fact that the vertical ozone and 
aerosol profile is generally not Known and the vertical 
water vapor profile would have to be interpolated from 
radiosonde observations at 12-hour intervals, and due to the 
circumstance that significant errors will probably occur 
only for very large zenith angles, no attempts have been 
made to account for that. An example taking account of _ the 
vertical profile of the aerosol mixing ratio can be found in 


a study by Dyer and Hicks (1968). 


3.2.5 Ozone transmissivity 

The transmissivity of solar radiation through the ozone 
layer of the atmosphere was calculated by Lacis and _ Hansen 
(1974) as 


uv 
2) 


T, =1-a, (um) -a, (um) (3.18) 


with 


a’(um) =0.02118um/(1+0.042um+0.000323(um) ) (3.19) 
a 0.205 
a,(um)=1.082um/(1+138.6um) = + 

0. 0658um/(1+(103.6um) ) (3.20) 


u is the thickness of the vertical ozone column, measured in 
cm at normal temperature and pressure and m is_' the optical 


air mass. 
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The error in these formulas is supposed to exceed 0.5% 
for um>10cm in equation 3.19 and for um>icm in equation 
3.20. With u being usually about 0.5cm the limits for a 0.5% 
error are m=20 (or: 9288°, according to equation 3.16) for 
a’'‘and m#2 (or: 0260°) for a’. Thus, the 0.5% error limit is 
practically never exceeded for the visible absorption term 
but is almost always exceeded for the ultraviolet 
absorption. 


The error in Ca incurred by the inaccuracy of equation 


3.19 can be estimated from equations 3.6 and 3.18 by 


pre) Cn ei IR 
Ada Dany OF 0 Tees re Coineay, 
For a typical arctic situation, 6=70° , u=0.5cm, and _ the 


precipitable water w=0.3cm. Therefore, m=2.9, um=1.45cm, 
a’’*=0.029, a. =0.022. As will be shown later, Teg =0.791 
and a,=0.095. Thus, for a 5% error in a,, the error in 7, is 
+0.0005, and for a 50% error in a, , the error in ‘a is 
+0.005. In other words, as a 50% error is not likely, the 
inaccuracy in the ultraviolet absorption term for um>1icm is 
most probably not capable of seriously affecting the t.- 
results. 

Ozone is measured once daily at Resolute by different 
methods (Canada, Environment Canada etc., 1975ff.). No 
interpolation scheme for intermediate hours was introduced 


in the computer program. Day-to-day differences in ozone are 


usually not greater than 0.03cm NTP. For u=0.5cm and m=2.9, 
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the error in ozone transmissivity is Al, =20.002. The error 


in @, can again be estimated by 


z nf lg Tr 
Qa dT. (ie m Terns ie (RSF oe.) 
For the conditions mentioned above, the error in %, would 


amount to about +0.0008 and is completely negligible. Now it 
is also clear that errors in the measurement of u_= and 
interpolated u-values for days with missing readings do not 


pose a threat to the precision of the results. 


3.2.6 Rayleigh transmissivity 

The Rayleigh transmissivity T, of the atmosphere has 
been calculated various times with increasing accuracy since 
Rayleigh discovered a theoretical formulation for the amount 
of light scattered by molecules in 1871. A recent attempt 
was made by Davies and Idso in 1979. The result is expressed 
as a function of optical air mass m, for optical air mass 
increments ofA m=0.2. 

For greater accuracy, linear interpolation was used in 
this thesis. Some Tp -values may show that the transmissivity 


decreases only slowly with increasing air mass: 
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multiplied by a pressure correction factor of p(kKPa)/101.3. 
The publication of Rayleigh transmissivities with 4 

significant digits and the fact that linear interpolation 

produces errors of well below +0.001 implies that the error 


in ©, stemming from T, is negligible. 


3.2.7 Water vapor absorptivity 

The absorptivity of solar radiation in air as a_ result 
of the presence of water vapor was calculated by Yamamoto 
(1962) from laboratory measurements made by Howard et al. 
(1955) (cited in Yamamoto, 1962). Lacis and Hansen (1974) 
found an empirical formula for Yamamoto’s curve’ which 


applies to the atmosphere: 
0.635 
ay=2.9wm/ ((1+141.5wm) +5 .925wm) 23) 


w is the precipitable water measured in centimeters and m is 
the optical air mass. 

As the measurements were made at standard temperature 
and pressure, and water vapor absorption is temperature and 
pressure dependent, a correction has to be applied (Lacis 


and Hansen, 1974): 
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troposphere, the choice of appropriate values for pressure, 
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p, and temperature, Te is a problem. Furthermore, the 
exponent, n, is only Known to be approximately between 0.5 
and 1.0. Lacis and Hansen, for example, take arbitrarily 
p=81kKPa and make their calculations for n=0, 0.5, and 1. 

As the uncertainties with these corrections are 
obviously great, no pressure correction was applied here 
and, as a first guess, the surface temperature was taken in 
the temperature correction term. 

For the same situation used above (0=70°,u=0.5cm, 
w=0.3cm) the error in a, due to the neglect of the pressure 
correction can be estimated using p=81kKPa and n=0.75: ay, 
turns out to be in error by aa, =-0.004. The influence of 


errors in a, on T, may be estimated by 


ee ae (3.25) 


which yields in this example AT,=0.002. As will be seen 
later, this error is still about one order of magnitude less 
than the errors’ introduced by the uncertainties of the 
radiation measurements, and shall therefore be neglected. 
Another error may be introduced by the use of equation 
3.23 itself whose uncertainty supposedly exceeds 1% for 
uw >10cm. The highest values of w occur in summer and may 
reach 1.5cm in the Arctic. With a high optical air mass of, 
e.g., m=10 (@=85°), wm becomes 15cm. u is typically 0.35cm 
in summer. Thus, T,=0.910, T,=0.611, and a,=0.199. FOR* sa 
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error ina,, the error in Gis +0.006 and approximately 
DM 7 abracton Of pie gGreaten athan pathe: yerrnor waincurred: : by 
neglecting the pressure correction. As the events of high 
precipitable water content and low solar elevation do _ not 
often occur simultaneously, the unknown error by 
extrapolating uw beyond the 10cm limit may not appreciably 
influence the final T,.-values. 

Precipitable water is being calculated twice daily for 
Resolute (for about 0600 hours’ and 1800 hours LAT) from 
radiosonde ascents. Linear interpolation is performed in the 
computer program fe) estimate precipitable water for 
intermediate hours. But precipitable water can change 
rapidly rather than smoothly with the advection of a 
different air mass (not to be confused with optical air 
mass). In a typical summer case, m=2.9 (or: 6=70°), w=1.0cm, 
and aw=0.3cm over 12 hours. The typical maximum error. by 
interpolation is therefore 0.3cm, and the maximum error’. in 
a, is then about +0.011. Using equation 3.25 with u=0.35cm 
the maximum error in @, amounts to +0.006. But the fact that 
abrupt changes in precipitable water do not occur too often 
makes this error negligible for T,-mean values. 

Precipitable water readings are presented with an 
accuracy of 0.01in (or: 0.025cm). This suggests that the 
errors from the published values are one order of magnitude 
less than the errors from the sudden jumps of precipitable 


water. 
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In some instances, precipitable water readings were not 
available. For these hours precipitable water was estimated 
from dewpoint temperature readings using a regression of the 


type 
In(w)=a+bt, (3.26) 


where t, denotes dewpoint temperature and 4@ and 6 are 
constants. This formula was first derived by Reitan (1963) 
for monthly values of w and t, and was_ shown to be 
applicable also for hourly values (Bolsenga, 1965). 


For Resolute the equation was found to be 


In(w) =-0.0878+0.0538t, oma it 


with w incm and ty, in °C. To establish equation 3.27, 182 
cases of hours with uninterrupted sunshine and simultaneous 
measurements ot w and ty, were used. The correlation 


coefficient: 7s 0793. 


3.2.8 Selection of sunshine hours 
The shortest periods for which climatic data are 
collected are hourly periods. Thus, to evaluate equation 
3.6, hours with uninterrupted sunshine had to be found. 
Sunshine observations are made with the Campbel1-Stokes 
sunshine recorder (Canada, Environment Canada etc., 1978ff. ) 


which is basically a glass sphere with a diameter of 10cm 
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burning a narrow paper strip. The solar radiation. flux 
density perpendicular to the rays must exceed values of 
0.2cal/(cm?min) (=140Wm’) (Bider, 1959) to 0.3cal/(cm*min) 
(=210Wm~*) (Brooks and Brooks, 1947) to give a burn. The 
threshold value is dependent on moisture and temperature and 
therefore tends to be higher for sunrise than for sunset 
(Bider, 1959). Taking, for example, %.=0.085, u=0.45cm, and 
~ w20.3cem, 140Wm~* and 210Wm *correspond to solar zenith angles 
of 86.6° and 85.3”, respectively. 

There are two possible errors associated with 
Campbel1-Stokes sunshine records: Firstly, the glass sphere 
may give a burn with the presence of clouds. To avoid this, 
all occurrences of uninterrupted sunshine observed in_ the 
presence of As, Cc, Cs, Ci, or fog at any cloud level were 
ignored. Secondly, if the solar disk is obstructed by clouds 
for short periods of about one minute or less the paper 
strip will show uninterrupted sunshine. Further errors are 
introduced by the fact that sunshine duration is published 
in tenths of an hour’ so that even a three minute 
interruption of direct sunshine would be recorded as a 
sunshine duration of 10/10. No remedy for these errors is 
available. 

It is obvious. that these errors systematically 
contribute to raising the t,-values. Another error which 
cannot be accounted for, but with the same tendency of 
raising the t,-values, stems from the pyranometers. Their 


readings are correct only if the glass hemispheres’ are 
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completely clear. Any precipitation, dew, or hoar_ frost 
systematically lowers the radiation income and increases the 
resultant @,-values. Therefore, unexpectedly high T, -values 
have to be dealt with great caution, and may not represent 
the quantity they are supposed to stand for: aerosol optical 
depth. 

For the three-year period from January ist, 1978 to 
December 31st, 1980, 3183 hours of uninterrupted sunshine 
could be found at Resolute. But only for 537 hours are all 
necessary data available and no clouds of the types As, Cc, 
Cs, Ci, or fog reported. Of these 537 hours, 313 hours were 
reported to have 0/10 cloudiness, the remaining 224 hours 


had partial cloudiness of the cloud types not excluded. 


3.3 Results and preliminary discussion 

Appendix A contains all 537 hourly input data as well 
aisethe catculations Of cosG; 1T,, I, j1as, and=tee% eFigure 3 
presents the monthly results of 537 calculations of the 
aerosol optical depth at Resolute. As the distribution of 
the T,-values is usually not symmetric, medians were plotted 
rather than arithmetic means to avoid undue weighting of 
some, maybe erroneously, high ‘t.-values. Section 3.2 has 
shown in detail that only two errors are important in_ the 
calculation of tT, from equation 3.6: The errors in the 
pyranometer readings, which were reported as being as_ high 


as 5%, and the errors in the solar constant, being 
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approximately 14%. 

To account for these errors all calculations were 
performed again with all pyranometer readings lowered by 5% 
and the solar constant simultaneously raised by 1%. Then, 
al] calculations were performed with all pyranometer 
readings raised by 5% and the solar constant simultaneously 
lowered by 1%. The first procedure results in higher, the 
second in lower t.-values. These adjusted medians are also 
shown in Figure 3. 

Each error bar in Figure 3 contains a wide and a narrow 
section: The wide section indicates the range between the 
median of ©, with all pyranometer readings lowered by 5% and 
the solar constant raised by 1%, and the median of T, with 
all pyranometer readings raised by 5% and the solar constant 
lowered by 1%. The narrow section indicates the range within 
which 68% of all 7,-values can be found jin the original 
calculation with the original pyranometer readings and a 
solar constant of 1353Wm~. 

In Figure 3 only those months are shown for which at 
least 10 calculations of ta could be made. Months with fewer 
than 10 calculations were not rejected if data were 
available for at least 3 different days. 

The monthly results averaged over three years can be 
expressed as follows: 

Feb. 0.04 
Mar. 0.08 
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The results show a systematic seasonal march of the aerosol 
optical depth: ©, increases from low values in late winter 
to a distinct maximum in April and May. It decreases in June 
and remains at a medium to high level all summer. 

The very high values for August 1980 are probably due 
to the eruption of the volcano Mt. St. Helens in Washington 
State on May 18th, 1980. Similar repercussions of volcanic 
events on the aerosol optical depth have been studied many 
times. An impressive example for the diffusion of high 
aerosol turbidity after the eruption of Mt. Agung on Bali on 
March 17th, 1963, was documented by Dyer and Hicks (1968) 
and Volz (1970). In this case the aerosol optical depth 
reached its normal value only after more than two years. 

To elucidate the volcanic influence on the aerosol 
optical depth at Resolute, the medians for June, July, and 
October 1980 were plotted with dots in Figure 3. They were 
not plotted with solid lines because they do not fulfill the 
criteria for medians set above. The volcanic influence 
reached its maximum in August 1980 with a delay of 3 months 
and tapered off in fall 1980. But this is only true if the 
seasonal ‘background aerosol’ itself (solid line in Figure 
3) remains constant in fall which cannot be inferred’ from 


Figure 3. 
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Setting aside the unique volcanic event in 1980, the 
following can be said about the aerosol optical depth at 
Resolute compared to other arctic stations (see Tables 2 and 
ot 

CT. is approximately of the same size as in other arctic 
locations. This is especially true for summer; the values 
for the period of highest turbidity, spring, are somewhat 
lower at Resolute. This is one more confirmation that the C, 
-values calculated from integrating nephe lometer 
measurements of [as (bracketed in Tables 2 and 3) are much 
too low. 

The results from Resolute also confirm that it is not 
justified to call winter the season of high and summer the 
season of low turbidity. At least for Resolute, spring and 
only spring is the season of high turbidity whereas winter, 
as far as can be seen from Figure 3, is rather a period of 
low turbidity. The arguments in favor of or against a 
particular continent as the origin of high turbidity in the 
Arctic must be seen in the light of these findings. But a 
detailed discussion of the problem of the origin must be put 
off until further computational results are Known in Chapter 
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3.4 Further results 

It might be of interest to find out whether 7, for 
cloudy hours is different from 1, for cloudless hours. For 
cloudy hours t, might be greater, because with the presence 
of clouds there might be some atmospheric layers, even in 
the cloudless portion of the sky, where relative humidity 
approaches 100%. 

In recent years, a multitude of studies have been 
published (e.g. Covert et al., 1972; Hamel, 1972; Zuyev et 
al., 1973; Hanel and Bullrich, 1978; Fitzgerald, 1980) which 
examine the functional relationship between the aerosol 
particle radius (or the scattering coefficient) and relative 
humidity. As a rule it can be said that the particle radius, 
and hence the scattering coefficient, increases only 
slightly as long as relative humidity is below about 804%. 
Above about 80%, particles grow very fast and reach the size 
of small cloud droplets at 100%. This is the well-Known 
phenomenon of haze. 

The median of %, for all 537 cases considered is found 
to be 0.088. For the 313 cases of cloudless skies Ta is only 
slightly different and amounts. to 08 087% Also, the 


ic 


correlation coefficient between ©, and cloud amount is found 
to be only 0.26. Thus, the effect of increasing Tt, with 
increasing cloudiness could not be proved. 

The correlation coefficient between T, and relative 
humidity at station height amounts to only 0.06 and also 


fails to show the desired relationship. This is so because 
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relative humidity exceeded 80% only in 17% of all cases. 
Visibility is being reported hourly at Resolute. It 


co 


should be expected that vt, is high for low visibilities and 
vice versa. But the correlation coefficient amounts to 0.01 
and is even of the wrong sign. This is probably due to _ the 
fact that visibility is mot measured, but estimated from 
watching some remote objects (mountains etc.) at Known 
distances. ieehen aor ee visibility is not a continuous 
variable as only a limited number of landmarks are used _ for 
its estimation. 

An important question is whether @, shows’ different 
values for winds from certain directions or for calms: 

To examine the dependence on wind direction winds’ were 
divided into two classes: northerly and easterly winds 
blowing from directions S45. Are clockwise through 135°, and 
southerly and westerly winds from directions between 135° 
and 315°. To determine whether occurrences of very high 
aerosol turbidity (T,>0.168) or very low aerosol turbidity 
(t,< 0.056) are associated with winds from a= particular 
direction X-tests were performed on Tables 7a and 7b. The 
results were negative, or, in other words: No- significant 
relationship could be established between occurrences of 
high turbidity and wind direction, nor could a_e significant 
relationship be established between low turbidity and winds 
from a particular direction. 

For Resolute, the northerly and easterly directions are 


those facing Asia and Europe, and the southerly and westerly 
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directions are facing North America. The lack of any 
relationship is in stark contrast to the claims of Peterson 
et al. (1980) that northerly winds at Barrow carry higher 
aerosol turbidities than southerly winds (see Table 2). 

To examine whether there is any relationship between 
the occurrence of calms at Resolute and the occurrence of 
high turbidity (t,>0.168) or the occurrence of low turbidity 
(T, << 0.056) two more X -tests were performed (Tables 7c and 
7d). Once again, the results were negative indicating that 
occurrences of extremely high or low turbidity are not 
linked to calms. 

Only the median of the t,-values at calms is slightly 
higher than the median for all cases: 0.103 versus 0.088. 

From this it must be concluded that high turbidity is 
not a local phenomenon, i.e., not induced by local pollution 
sources. But a more thorough discussion of these findings 
has to wait until, in a later chapter, more computational 


results are Known. 


he 9 at 


: ; vw, = ni * it 7 


eae a oh 


) (att ore a 4s enot 
yrs vos ‘ ai hoigstis sae 


noenstet Fo embets orth of Sesrsnas 
seriptd qaiss ened ote ina ‘ +4 
1S efdsT eae) aborntw inten’ v 
neewied didenctistey yns st sre sata artnet 
to soneraNosd oft bre ehitoae mite anit go! Foo gare utso 
e+ tbidwi wot 7 SIS WGIS ert 46. {8b 20, a) wtibtdawd a 
bn ov gehdisT) bento tte | aibw eiael~ be sic Ow) \s20098 
fer? gntteotont evi taper’ orew at 1ded oct Atos e000 8 
ton ae. yitbtd«wt wool 3% — viomend xe +g a@eonet4 
| beg ee vembe ot wolat 

Viinpilea- et elise te iemitente ant to netbem ent ube hiar | bigs 
BHO.0 euesev 601.0 resesd Peto? netbem ott hart seripid 

2 yirotdews dota tert bat grind Sc teum the etrty’ an 
notiviteg toot xd besubnt farnr),.2.' nonsmunariq (soot ¢-ton — 


“ae 
soritont sees Yo not eayoe tb niques encom « tue a) vat 
‘enciaieans ecm retqsria aoet s mf fron ttew oF . a 

~ eee 846 estueer 


a : em & 


CO16S > O05 


(c) 


es < 0.056 »0.056 


(b) 
winds 89 386 


calms 


(d) 


58 


Table 7: Simultaneous occurrences of high (a) and low (b) aerosol turbidities and winds 
from either northerly and easterly or southerly and westerly directions. Simultaneous 
occurrences of high (c) and low (d) turbidities and calms. 
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4. Aerosol single-scattering albedo at Resolute 


4.1 Introduction 

There are var 10uUSs methods to determine the 
single-scattering albedo or, alternatively, the imaginary 
part of the refractive index of an aerosol. They can be 
divided into two classes: laboratory methods and 
field-experiment methods. 

The majority of the results given in Tables 4 and 5 
were obtained in laboratories using atmospheric aerosol 
samples collected either at ground or by aarcratt. 
Laboratory methods, in turn, can be divided into’ two 
classes: Either the volume-absorption and volume-scattering 
coefficients are determined directly or the angular 
distribution of the scattered radiation is recorded with a 
scanning nephe lometer and the result compared with 
scattering phase functions derived from complete solutions 
OF the Mie scattering theory for various particle 
size-distributions and imaginary parts. 

There are three major disadvantages associated with 
laboratory methods: Firstly, there is reason to fear that 
the aerosol undergoes changes , e.g. in size-distribution 
and water coating, after it has been collected. Secondly, 
samples taken at ground level do not represent the mean 
conditions of the entire atmosphere, only those of _ the 


surface boundary layer. Thirdly, since laboratory methods 
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are time and money consuming, the number of samples is 
generally not sufficient to gain climatological mean values. 

Field-experiment methods, on the other hand, directly 
measure the interaction of solar radiation with the aerosol 
in the atmosphere. Therefore, they are in essence’ solutions 
to the problem of solar radiation transfer. As the method 
used in this thesis can also be classified as a 
field-experiment method, a short literature review of other 
similar approaches will be given on the following pages. 

Roach (1961) and Method and Carlson (1982) made 
radiation measurements from aircraft at different altitudes 
and immediately derived the absorbed, upward scattered, and 
downward scattered fluxes in an atmospheric layer.Generally, 
the use of aircraft immediately provides the data necessary 
to calculate w,. But, for obvious reasons, #4 has to be 
determined from ground-based measurements in most cases. 
This implies that models have to be used that predict’ the 
diffuse downwelling radiation at ground level. 

Forbes and Hamilton (1971), used a very simple 
radiative transfer model which gives doubtful results, as 
noted in Chapter 2. 

Herman et al. (1975) present the so-called 
diffuse-direct radiation method: It is based on_= spectral 
intervals for which only Rayleigh scattering and aerosol 
extinction are relevant. The radiative transfer problem is 
solved using a plane-parallel atmosphere with a presumed 


vertical aerosol distribution according to Elterman (1964) 
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and a Junge aerosol size-distribution. Calculations of the 
total diffuse flux at the ground Dl as a function of the 
solar zenith angle 0, the wavelength A, the aerosol optical 


depth t,, the complex index of refraction n=n,-in,, and the 


/ / 


surface albedo ‘a’ are thereby obtained. The direct’ solar 
flux at the ground on a surface perpendicular to the solar 
rays I” can also be obtained (see Chapter 3). Next, a 
multitude of diagrams is plotted showing the ratio DY / 1’ 
versus n, with 9, A, T,, and ‘a’ as parameters. With DL and 
1° measured for a wavelength A, and 9, T,., and ‘a’ known, n, 
can be derived by interpolating between the diagrams. This 
method has also been used by DeLuisi et al. (1976) and, ina 
modified version, by Carlson and Caverly (1977). The 


foregoing diffuse-direct method was extended by King and 


/ / 


Herman (1979) to derive n, and ‘a’ simultaneously and was 
applied by King (1979). 

The diffuse-direct method cannot apply to the data 
available for the Canadian Arctic because spectral radiation 
measurements are not taken there. Another disadvantage of 
the 1975-version is the fact that interpolations have to be 
performed between various diagrams. The 1979-version, 
according to King (1979), necessitates ‘a clear atmosphere 
for a long enough period of time for the solar zenith angle 
to undergo a_ large change, preferably of about AC ee alias 
requirement also excludes its application to the Arctic. 


The foregoing introduction has introduced and 


eliminated approaches which calculate the single-scattering 
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albedo using airborne measurements or spectral measurements 
because such measurements are not available for Canada’s 
APG Lic. 

The present author believes that no attempts have so 
far. been made to calculate w, or ne in the most 
straightforward way which is to solve empirical or 
theoretical radiative transfer models for © or n,. 

In this study the present author attempts to apply 
increasingly complicated models for all hours’ with a 
sunshine duration of 10/10 and 0/10 cloudiness (and all 
other relevant data available) to calculate w, for’ the 
entire solar spectrum at Resolute for the period from 


January 1978 to December 1980. 


4.2 Davies-Hay model 

There have always been two largely independent 
approaches to the radiative transfer problem: Firstly, there 
is a purely empirical approach which might also be called 
the single-scattering approach, because it assumes that 
diffuse radiation is only scattered once on its way from the 
sun to the bottom of the atmosphere. To use this approach 
transmissivities for the whole atmosphere are needed and, 
for the most part, readily obtained. Secondly, there are 
theoretical or multiple-scattering approaches, a very 
popular version of which is the two-stream approximation. 


For this approach, volume-extinction coefficients are needed 
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which are, for the case of the entire solar spectrum, not as 
easily obtained as are the transmissivities used in_ the 
single-scattering approach. 

Sheppard (1958) suggests that, as a thumb~ rule, 
single-scattering approaches may be used if the _ total 
transmissivity of the direct solar beam exceeds 0.90, and 
multiple-scattering has to be used if it is less than 0.74. 
This suggests that in the Arctic, with its low Rayleigh 
transmissivity, the latter is always mandatory. 

This section deals with the latest version of the 
single-scattering approach, that of Davies and Hay (1980), 
whereas the following sections deal with various’ two-stream 
approximations. 

The grandfathers of empirical radiative transfer are 
Linke and Angstrom, for total and spectral direct solar 
radiation fluxes, respectively. Linke’s turbidity factor was 
used from the 1920s to the 1950s. In 1954, Houghton came up 
with a new empirical radiative transfer model which was 
still rather simple with respect to the atmospheric aerosol. 
This model was modified by Monteith (1962), Idso 
(1969,1970), Davies et al. (1975), Suckling and Hay (1976), 
Paltridge and Platt (1976,p.121), Davies and Idso (1979), 
and finally by Davies and Hay (1980). This latest version is 
a sophisticated empirical solar radiation transfer 
approximation, and will be used as ae first attempt to 
calculate the single-scattering albedo. Hereinafter it will 


be referred to as Davies-Hay model. 
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According to the Davies-Hay model the downwelling 
diffuse radiation Dl) at the bottom of the atmosphere for 


cloudless conditions can be written as: 


D) =D, +D, +D. (4.1) 
with: 

D, =7, —+* T, 1, cos@ (4.2) 

Dwerilpatercays eet Jw, Bg 1, cos O (4533) 

D, =ax(1+Dg+D,)/(1-ax) (4.4) 


De ,D, , and Ds are the diffuse fluxes due to Rayleigh 
scattering, aerosol scattering, and multiple reflection of 
global radiation between the surface and the clear 
atmosphere. “respectively. “lay ry °O, 1), “To 1 "and ay ‘have 
the same meanings as_ in the previous chapters. 4 is the 
single-scattering albedo for the aerosol, B, is the ratio of 
downward to total aerosol scatter for direct radiation, a is 
the surface albedo and « is the albedo of the clear sky for 


upwelling global radiation which is in the Davies-Hay model 
=0.0685+(1-T, )w,(1-B! ) (4.5) 


The first term represents the albedo due to Rayleigh 
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backscattering and the second term the albedo for aerosol 
backscattering. Bt is the ratio of upward to total aerosol 
scatter for upwelling global radiation. 

The surface albedo can be obtained from measurements of 
the upwelling global and the downwelling global radiation. 
Both are measured at Resolute. The only new quantities are 
B, and Bl. 

B. 1s a’ funetion of the’zenith angle of the “impinging 
radiation. This is so because the aerosol scattering phase 
function (which indicates the amount of radiation scattered 
into an infinitesimal solid angle as_ function of the 
scattering angle) is highly asymmetric with respect to the 
scattering angle but symmetric with respect to the azimuthal 
angle and has the shape of a club with a marked forward 
peak, minima at around right angles to the impinging 
radiation, and a minor maximum in the backward direction. 
Now, it can easily be understood that with the sun being at 
zenith, ‘the club points downward’, in other words, most of 
the radiation is scattered downward to the earth. At sunrise 
and sunset, ‘the club is lying horizontally’ and, therefore, 
whatever shape the club has, 504 of the radiation is 
scattered down to the earth and 50% back to space. 

Various values for the percentage of downward to total 
aerosol scatter, Ba, with the sun being at zenith have’ been 
measured: 0.94 (Roach, 1961), 0.92 (Robinson, 1962), 
0.80-0.94 (Charlson et al., 1974), 0.83-0.94 (Waggoner et 
al., 1981), and 0.88 (Method and Carlson, 1982). Much more 
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Gittreu)t] to rind in literature, is, the, variation. of Bawwath 
the solar zenith angle 9. So far, only Robinson (1962) has 


published measurements: 


OR DS On 0.6 O55 0.4 O13 


Therefore, these data were taken as an estimate and 
graphically interpolated for all values of cos@. As B, (0=0" ) 
=0.92 seems to be slightly too high with regard to the more 
recent findings, B,(9=0°)=0.90 was taken and the whole 
function B,(8) proportionally shifted to slightly lower 
values. 

The upward to total aerosol scatter for upwelling 
radiation pf poses another problem: No uniform zenith angle 
can be allotted to diffuse radiation. But, as will be shown 
in the section on two-stream approximations, it is common 
practice to assume that, at least upwelling, diffuse 
radiation traverses a ‘'mean’ optical air mass of 1.66 which 
corresponds to a zenith angle of 53. Therefore, pt =B, (@=53° | 
=0.81 was taken. 

Three more things have to be explained in the 
Davies-Hay model (equations 4.1 to 4.5): 

Firstly, the scattering phase function for Rayleigh 
scattering has the shape of a bone with identical maxima in 
the forward and backward direction and a symmetrical minimum 


at right angles to the impinging radiation. Therefore, the 
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downward to total scattering ratio is always exactly 50% for 
any ‘location of the bone’. This is the reason for the 
factor 1/2 in equation 4.2. 

Secondly, as Rayleigh scattering is restricted to the 
shorter wavelengths and water vapor absorption to the longer 
wavelengths, no water vapor absorption term was’ introduced 
in equation 4.2. 

Thirdly, the multiple reflection term (equation 4.4) 
can easily be understood after a Taylor-expansion is applied 


to it. This yields: 
D.=(1+D,+D, )(ax ta ata w+, ..) (4.6) 


Now, equations 4.1 to 4.5 can be solved for w, because 
all other variables are Known, namely: 
Jit rOmeeduat1Onsio 1G: Lor, 20; 
Siem th On S6C MlOl (ont 203; 
det romeduatlOonicc.20. 
-T, =exp(-T,m): T, from the result of calculations made in 
Chapter 3; m from equation 3.16; 
=1,.: from™=seeriorn 372.41; 
-cos@: from equations 3.9 to 3.12; 
-B, BI: from Robinson (1962), as discussed in this section; 
-a : from measurements of the global downwelling and 
upwelling radiation. 

The results of these calculations will be given at’ the 


end of this Chapter. 
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4.3 Two-stream approximation 

The two-stream approximation of radiative transfer was 
proposed by Schuster (1905) and has since been subjected to 
considerable sophistication. The field of diffuse solar 
radiation at some height z in the atmosphere is thought of 
consisting of two ‘streams’ or fluxes: the upwelling and the 
downwelling diffuse flux, Dt(z) and Dl(z). In a _ height 
element dz, the two fluxes are subject to changes dDf(z) and 
dDlLiz), respectively. They decrease because of absorption 
and scattering-out, but they increase because part of the 
opposing flux is scattered in and because at any height of 
the atmosphere the direct solar beam acts as a source of 
diffuse radiation through scattering. 


This yields a system of two inhomogeneous differential 


equations: 
ATE) 2-4, Dt (z)+a,.Db (z) +b, f(z) 0477) 
APUG) --2, Dt (z)+a,D Uz) -b F(z) (4.8) 


If we assume that ozone absorption is already completed 
when the solar beam hits the top of the atmosphere, only 
Rayleigh scattering, aerosol scattering, water vapor 
absorption, and aerosol absorption have to be condidered. 
Then the matrix elements a;,, the vector elements b;, and the 


source function f(z), describing the direct beam, become: 
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a, emt (e+tpt+ A, (1-82) ) (4.9) 
a, amare wrbticB ile (4.10) 
a, =mtit fit wetllias De (4. v1) 
at rai ESP Palace Tale (4.12) 
Beal, ace One Be))) Raed 
be ily ipl a4ie* fe, B,) (4.14) 


Fiz. -alpld, = 


=(exp(-T,m)-(1-exp(-T,m) ) )exp(-T, m) = 


= (exp ( -3, 2, ) +exp (0, 452 )-1)exp(-A4et) (4.15) 


cos 9 


At this point it should be mentioned that in all 
previous two-stream approximations Known to this author the 
distinction between the various terms involving atmospheric 
absorbers and scatterers in equations 4.9 to 4.15 were not 
as thorough. The problem of evaluating these coefficients 
will be dealt with in more detail in Section 4.4 than can be 
found elsewhere in literature. It must therefore be stressed 
that the set of equations 4.9 to 4.15 is based on a physical 


development more so than on previous literature. As_ the 
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inhomogeneity or source function (equations 4.13 to 4.15) 
has not been used in this form before, a complete solution 
of this system of inhomogeneous differential equations could 
not be found elsewhere in literature. 

Let us now return to equations 4.9 to 4.15. mf and mJ] 
are the diffusivity factors for upwelling and downwelling 
diffuse radiation, respectively. Details of their meaning 
will be given in Section 4.4.1. 

Be ice By, and Ty are the water vapor vo lume 
absorption coefficients and the water vapor optical depth, 
respectively; Bi, py, Ba, and T, are the Rayleigh volume 
scattering coefficients and the Rayleigh optical depth, 
respectively; the definitions of these quantities are 
analogous to the definitions of a and T, (see equations 2.9 
and 2.4). Reasons to introduce three different Bw 3 
coefficients and three different /3g -coefficients will be 
given in Sections 4.4.4 and 4.4.3., respectively. 

A distinction between pl and BY has also been made, the 
first variable denoting the upward to total aerosol scatter 
ratio for upwelling diffuse radiation, and the latter one 
denoting the downward to total aerosol scatter ratio for 
downwe 11ing diffuse radiation. The reason for this 
distinction will also be given in Section 4.4.1. 

A finite height H of the atmosphere has been introduced 
to find an expression for the source function f(z) in 
regular height coordinates. The implications of this step 


will be discussed in Section 4.4.2. All other quantities in 
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equations 4.9 to 4.15 are the same as defined earlier. 

Two things have to be said about the source’ functions 
b, f(z) and b,f(z), which describe the conversion of direct 
solar radiation to diffuse radiation. As T, is not 
incorporated in f(z) the total transmissivity for the direct 
beam is 1,(Tpg-a,)T, rather than (T,Tz-a,)T, as used in 
Chapter 3. As T, is almost always greater than 0.90 and in 
many cases close to 0.95 the difference between these 
formulas is typically of the order of 1%. Furthermore, it 
should be mentioned that the cos@-term does not appear in 
the source function. To make this plausible, consider a 
volume element in the atmosphere and imagine, for the 
moment, that the direct solar radiation streaming through 
this volume element has not been attenuated at all before it 
reaches the volume element. Then it is immediately obvious 
that the amount of radiation scattered inside the volume 
element, in other words the intensity of the source of 
diffuse radiation, is the same for any zenith angle of the 
impinging solar radiation. The fact that, in reality, 
radiation has already been attenuated before it reaches the 
volume element and that the intensity of the source is'- thus 
a function of the solar zenith angle is, of course, being 
considered by means of the factors 1,, T,, T., and a,. 

The analytical solution of this system of two 
inhomogeneous differential equations (equations 4.7 to 4.8) 
is tedious but mathematically straightforward (see, e.g., 


Bronstein and Semendjajew, 1976, pp.391f.): The most general 
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solution of the homogeneous system (i.e., for f(z)=0) turns 


out to be 
Dit zleGyexolr, 2) +c, expr, z) (74'14.6)) 
Bitzi — Chexpir)z)+22== Cl expir, <) (4; aii} 


where r, andr, are the eigenvalues of the homogeneous 


system 

rn, =(a,.7a, +) (a,,-a,,) -4(a,a,,-a, ar.) ) (4.18) 
Wwieen) theron IC OsandeCiy are ianbiitradr y-cons tan tsi-velo tgets one 
particular solution of the inhomogeneous system (i.e., for 


f(z)#0) we assume C,=C,(z) and C,=C,(z) and insert equations 
4.16 and 4.17 into 4.7 and 4.8. This yields an inhomogeneous 
system of algebraic equations with the unknowns G-(z) and 

/ 
C 


,(Z) where the prime denotes the first derivative with 


respect to z. This system can immediately be solved for 
Eira) and Cy(z), and C,(z) and C,(z) can be found by simple 
integration. As we need only one particular solution we set 
the integration constants equal to zero. Inserting C,(z) and 
Co Ange Teo cctiatrOncme 4 slO-sand: et) yields the desdmed 
particular solution. The most general solution of the 
inhomogeneous system is then’ simply the sum of the most 
general solution of the homogeneous system (equations 4.16 


and 4.17) and the particular solution of the inhomogeneous 
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system. 


This system still has two arbitrary. constants, C, and 


Co . Therefore’ two boundary: conditions "are required. “They 
are: 
DI Z2sh)=0 (4.19) 
DiahzsO.) san dez= 0) 41 z= 0)))) (oe 0) 


I(z=0) means the direct horizontal solar radiation at the 


/ 


ground and ‘a’ the surface albedo. 
As we are only interested in DJ(z=0) (in the Davies-Hay 


model denoted by DJ) only this result shall be given here: 


et NCR lee Ome nIe eT, (4.21) 
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Bela Diy cae Lie Due OE ea 


Qin Qin 
+cos OW exp(-F, H)+exp(-F, H)-exp(-F,H) (48 27h} 
F, =(Pe+ 8.) /cos @ (4.28) 
F,=(fu + ,)/cos @ (4.29) 
F, = 3./cos G (4.30) 


As D\Y is being measured, equation 4.21 can be 
comprehended as an extremely complicated transcendental! 
function in %, which can be numerically solved for w,. This 
requires that all the other variables occuring in equations 
4.9 to 4.15 are Known. This is actually the case. The next 


section will show how they can be evaluated. 


4.4 Evaluation of the two-stream approximation 


4.4.1 Diffusivity factors m and aerosol scattering ratios B 
The diffusivity factors, mt and ml, for diffuse 
radiation are ana logous to the optical air mass m for direct 
radiation. They are simply magnification factors, indicating 
by how many times the path length of radiation through the 
entire atmosphere is greater than the shortest possible path 


length at vertical incidence. For direct radiation this 
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problem is easily solved (equations 3.15 and 3.16). For 
diffuse radiation, whether upwelling or downwelling, = an 
appropriate magnification factor has to be found which jis an 
average value for radiation from one entire hemisphere. 

If the diffuse radiation is isotropic (i.e., the 
intensity is independent of direction for the entire 
hemisphere considered) then the problem can be_= solved 
numerically (see Kondratyev, 1969, p18f.) with the result 
that the diffusivity factor ranges between 1.2 for an 
intransparent atmosphere and 2.0 for a comp letely 
transparent atmosphere. But for the wide range of diffuse 
radiation transmissivities from 0.2 to 0.8 an average value 
of 1.66 can be applied with no more than 10% error. The 
factor 1.66 has found its way into many radiation studies. 

The problem of finding a diffusivity factor for 
isotropic diffuse radiation has also been tackled by 
Chandrasekhar (1950,p.54f.) who comes up with the result \3= 
=1.73... as a first order approximation to the problem. 

In recent years, the latter value has become popular 
(e.g., Sagan and Pollack, 1967; Liou, 1974; Liou, 1980, 
p.184ffis> Paltridge and Platt, 1976, p.100), although in 
some cases even the value 2.0 finds application (Paltridge 
and Platt, 1976, p.73f.; Chylek and Coakley, 1974). 

Unfortunately, the problem of finding appropriate 
diffusivity factors is much more complicated because the 
isotropy assumption does usually not hold. This is 


especially true for downwelling diffuse radiation. The 
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reason for anisotropy of downwelling diffuse radiation is 
the deviation from sphericity of the Rayleigh and, above 
all, the aerosol phase functions which were previously 
described as a bone and a club, respectively. Therefore, the 
bulk of downwelling diffuse radiation has a zenith angle 
similar to that of the sun. For upwelling diffuse radiation 
at height z, the isotropy assumption is fulfilled if the 
ground acts as a perfectly isotropic rather than partially 
specular reflector and the anisotropic upwelling diffuse 
radiation stemming from backscattering inside the atmosphere 
between the ground and height z is ignored for a moment. But 
it is Known from everyday experience that certain surfaces, 
above all flat water but also ice and melting snow covered 
with a thin ice sheet, give highly anisotropic and almost 
specular reflections. 

For some unkKnown reason the problem of anisotropy of 
diffuse radiation has largely been ignored in the discussion 
of diffusivity factors. Only  Kondratyev (1969,p.207) 
mentions some work done by Kondratyev'= and Sender ikhina 
(1959) to evaluate the diffusivity factors from intensity 
measurements over the upper and lower hemisphere. 
Measurements were made in the Leningrad region and on_ the 
Crimea but the surface conditions were not reported. Their 
result, formulated below, shows a _ considerable functional 
dependence of the diffusivity factors m1 and ml on the 
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mt=(0.21/cosQ)+1.73 (0<@<75° ) (A. 24) 
mv=(0.54/cos@)+1.09 (0<9<73° ) (4.32) 


For great zenith angles ch, which are of particular 
importance in the Arctic, these diffusivity factors may 
assume very great values and are considerably different from 
the diffusivity factors mentioned before. Therefore, it 
appears necessary that all calculations in this’ study be 
made with at least two different sets of diffusivity 
factors: In a first calculation the diffusivity factors of 
equations 4.31 and 4.32 were taken. These calculations wil] 
hereinafter be referred to as Kondratyev-model. In a second 


calculation the diffusivity factors 
ml=eml=1.66 (4.33) 


were chosen which will be referred to as 1.66-model. 

Now the appropriate choice of the upward to _ total 
aerosol scatter ratio pi for upwelling diffuse radiation and 
the downward to total aerosol scatter ratio BY for 
downwelling diffuse radiation can be made by assuming that 
the diffuse radiations from the lower and _ the upper 
hemispheres behave as_ if they had zenith angles Ot and OV 
defined by 
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@)=arccos(1/my\) (4.35) 


Then the aerosol scatter ratios Bl and BY can immediately be 
taken from the table given by Robinson (1962) (see Section 


4.2). 


4.4.2 Height of the atmosphere and volume = extinction 
coefficients 

The factors /3,, Ba» and Ba in equations 4.9 to 4.15 are 
not functions of height. Therefore, we have to introduce two 
assumptions: Firstly, we have to assume that the mixing 
ratios of water vapor and aerosol are constant for the 
entire atmosphere. In the next section (4.5) this somewhat 
unrealistic assumption will be relaxed by dividing the 
atmosphere into two layers at 50KPa and distributing water 
vapor and aerosol in a more realistic fashion. But even with 
the assumption of constant mixing ratios the vo lume 
extinction coefficients ws Be. and Ba taper off with 
height at the same rate density decreases with height. 
Therefore, an atmosphere with constant density has to be 
introduced which has a finite height. To relate this finite 
height with parameters of the real atmosphere an isotherma|] 
atmosphere has to be assumed (this assumption holds in the 
Arctic better than anywhere else). Integrating the 
density-height function of the rea | but isotherma | 


atmosphere 
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Q(z) =9(z=0)exp(-gz/RT) (4.36) 


with respect to z from 0 tow yields that an atmosphere of 


constant density 9=9(z=0) has a finite height H given by 
H=RT/g (4.37) 


where R is the gas constant for dry air, T the mean 
temperature and g the gravity acceleration. H is often 
called the scale height of the atmosphere. 


Now the volume extinction coefficients are simply: 


Pa=Te/H (4.39) 
fa =T./H (4.40) 


At this point two arguments might be put forward 
against the introduction of a constant density atmosphere: 

Firstly, it might be argued that the choice of any 
other density, e.g. 9=9(z=0), would have led to a _ height 
H/y. A brief check of the matrix and vector elements 
(equations 4.9 to 4.15) shows that the’ factor y would 
equally appear in all of them so the system of differential 
equations would essentially be the same except for the fact 


that any height z would now correspond to a height ZAK} - To 
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avoid this problem all heights have to be measured in’ terms 
of the finite height of the atmosphere H. This procedure is 
similar to the introduction of pressure coordinates. 

Secondly, it might be conjectured that the atmospheric 
radiation field would be different (in terms of the newly 
defined relative height coordinates, or in terms of pressure 
coordinates) for different density-height functions. This is 
certainly not true as all absorption and scattering 
processes are independent of each other and the total amount 
of absorption and. scattering depends only on the total 
amount of absorbers and scatterers in the atmosphere. 

The limitations of the independence of the amount of 
absorbed and scattered radiation on the density-height 
function are the following: For very high densities Rayleigh 
and Mie scattering theories are not applicable as they are 
based on the far’ field solution of dipole radiation. For 
very low densities the sphericity of the terrestrial 
atmosphere would bring about different radiation fields. 

To calculate H, an atmospheric mean temperature T has 
to be found for Resolute. From Palmén and Newton (1969,p.3) 
a density weighted mean temperature of the atmosphere at 
75°N can be calculated as being 233K in January and 253K in 
July. To interpolate for the intermediate months a_e simple 


cosine function was used giving the following results: 
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There is one more problem with the calculation of the volume 
extinction eceffictientse/3: Formulas 4.38 to 4.40 suggest 
that the coefficients [2 can directly be obtained from the 
respective optical depths © which, in turn, can be obtained 


from the respective transmissivities T, j.e., 

T=exp(-tm) =exp( -/Hm) (4.41) 
or solved for /3 

(B=-qeIn(T) (4,42) 


Here, m describes the optical air mass or the diffusivity 
factor, according to the type of radiation. 

But formula 4.41, commonly’ Known as Beer’s law, is, 
strictly speaking, only applicable for monochromatic 
radiation. As the direct beam transmissivities for the 
entire solar spectrum T of Rayleigh scattering and water 
vapor absorption are well known (see Section 3.2) we can 
insert some values into equation 4.42 and find that Br and 
[Sw are not constant but decrease with increasing optical air 
mass m. For H=7100m (corresponding to T=243K) and 2.3mm of 


precipitable water (a typical value at Resolute) we find: 
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The /3-values calculated with equation 4.42 are thus. only 
averages over the entire path length m. 

The reason for the /3-factors to decrease with 
increasing air mass is the following: Water vapor absorption 
is. shighily» selective, aia@+y.nl Eapoccurs:.only-,,intecentain 
spectral bands. Rayleigh scattering, although it is not 
selective, behaves in a similar way insofar as_~ shorter 
wavelengths are much more sub jected to it than longer ones. 
The result is similar: After radiation has largely’ been 
extinguished in the spectral regions preferred by the 
extinguishing processes only little is left te be 
attenuated. Therefore /3, which is according to equation 2.8 
the fraction of radiation absorbed or scattered per unit 
length, decreases drastically. 

As a transmission function T, for aerosol extinction, 
similar to the transmission functions Tp, for Rayleigh 
scattering and T, for water vapor absorption, could not be 


found, equation 4.40 was used to calculate [a 


4.4.3 Rayleigh volume scattering coefficient 

Equation 4.42 can be used to determine the Rayleigh 
volume scattering coefficient [3 for direct solar radiation 
in the source function (equations 4.13 to 4.15). But a 
considerable error would be introduced if the same were 
taken. in the elements of the matrix a; describing the 


diffuse radiation. fluxes. This -is due to the fact that 


diffuse radiation has a completely different spectral 
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composition than direct radiation. 

Kondratyev (1969,p.364) gives the theoretically 
calculated spectrum of downwelling diffuse radiation in a 
cloudless atmosphere. The mean Ray leigh scattering 


coefficient fx for this spectrum can then be calculated as 
BY=E(D1, PRIS, FEMA (4.43) 


where the summation goes over a finite mumber of = spectral 
intervals i, I; is the spectral flux in arbitrary units of 
the interval centered at A; and € is a constant which can be 
calculated from the information given by Paltridge and Platt 


(1976,p.119f) as being 
- 30 ~~ 3 
C= 054108 (273/ 1) (6/100) m (4,44) 
<1/ASturns out to be 4.36*10'm "so that (3x becomes 
(¥=4.56*10°°(273/T) (p/100)m”' (4.45) 


T and p are the surface temperature in K and the surface 
pressure in KPa, respectively. 

To calculate ay the Rayleigh volume scattering 
coefficient for upwelling diffuse radiation, we assume that 
upwelling diffuse radiation at any height z is exclusively 
composed of downwelling global radiation reflected at the 


ground and that the ground is a grey reflector. Now we have 
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to trace two radiation fluxes: the flux consisting of 
reflected direct radiation and the flux consisting of 
reflected diffuse radiation. The first of these two fluxes 
traverses the air mass m before it reaches the ground and 
the air mass m1 from the ground back to the top of the 
atmosphere. Therefore, its transmissivity on its way back up 


through the atmosphere is given by 
Ta (m+mt) /T, (m) i (4.46) 


and thus its respective volume scattering coefficient is 


given by (see equation 4.42): 


| Tr (mt mi) 
= —— |i 4.47 
where the superscript (dir) denotes its relation to 
reflected direct radiation. 

The second of these fluxes has a spectral composition 
which is identical to the spectral composition of 


downwelling diffuse radiation; therefore 
Gemeesee (4.48) 


where (dif) denotes reflected diffuse radiation. 
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4.4.4 Water vapor volume absorption coefficient 

Equation 4.42 can also be used to determine the water 
vapor volume absorption coefficient [Pw for direct solar 
radiation in the source function. 

To calculate os for downwelling diffuse radiation the 
assumption is made that the scattering processes’ which 
generate DJ (z=0) take place at z=H/2 so that it has to 
traverse an air mass of (m+ml)/2. Then equation 4.42 can be 
used giving 


i | 


Bre ayaa) (eT ea mem s)) (4.50) 


For the computation of fi), the water vapor volume absorption 
coefficient for upwelling diffuse radiation, the separation 


of the two fluxes made in Section 4.4.3 gives 
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tT (dir) ip (dif) 
aTeT Eze age (4.53) 
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Results of the Kondratyev-model and the 1.66-model will be 


given in the end of this chapter. 
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4.5 Extension to a two-layer atmosphere 

It has already been mentioned that the assumption of 
constant mixing ratios of aerosol and water vapor throughout 
the depth of the atmosphere is not justified, as the higher 
layers of the atmosphere don’ t usually contain much of these 
attenuating substances. A vertical profile of the aerosol 
concentration is not available on any regular basis but, as 
a first estimate, the atmosphere can be divided into’ two 
layers of equal mass, and_ the percentage of aerosol and 
water vapor in each of these layers can be estimated using 
commonly accepted values. 

The calculational procedures for a two-layer atmosphere 
are basically the same as for a one-layer atmosphere, the 
only difference being that all calculations have to be made 
separately for each layer. All variables pertaining to the 
upper part of the atmosphere shall hereinafter be marked 
with a prime. There are now four arbitrary constants, so the 
Same number of boundary conditions has to be found. The’ two 
new ones arise from the condition that the upwelling as well 
as the downwelling diffuse flux must be continuous at the 
interface of the two layers (see, e.g., Shettle and Weinman, 
1970; Weinman and Guetter, 1972; Liou, 1975). 


This yields the following set of equations: 
i 
D|(z=H) =0 (4.54) 


Bil ¢z=0 )-oaD jez 0+ HVZ=0 }, ) (4.55) 
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D |(z=H/2)=0 |(zsH/2) (4.56) 

DN z=H/2)=DNz=H/2) (4.57) 
Again we are only interested in Dj(z=0)=2DJ: 
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Biz " 


(4.58) 


D, and D, and all the other’ functions needed to 
calculate D, and D, are identical to equations 4.24, 4.25, 


4.28, 4.29, and 4.30. 
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I,=E, (from equation 4.27) (4.67) 
I,=—, (from equation 4.26, but all variables primed) 
(4.68) 
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Again, as DJ) is being measured, equation 4.58 is a very 
complicated transcendental function in Ww, and can only 


numerically be solved for it. 
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4.6 Evaluation of the two-layer model 


4.6.1 Vertical distribution of aerosol and water vapor in 
the atmosphere 

Let us for a moment return to an isothermal atmosphere 
whose density decreases in an exponential way, with H being 
the scale height of the atmosphere. Let us then assume that 
the attenuating constituent i (aerosol or water vapor) 
decreases also exponentially with height, but with a scale 
height He: Eet us ‘now®calculate the fraction t-R. vofis the 
total amount of the constituent i which can be found in the 
lower half of the atmosphere. 

If M; (z) denotes the mass of constituent i per height 
interval dz then the total mass of constituent i in the 
entire atmosphere is 


20 20 
{w, (z)az=m, (220) fexp(-2/H, dz=Hs Mi. (z=0) Ras 7y) 
i 0 

If our isothermal atmosphere has a temperature of 243K, then 

H=7100m and the height of the 50kKPa surface which separates 

the upper from the lower half of the atmosphere is found to 

be at 4900m. Now the total mass of the constituent i between 

the ground and 4900m can be calculated as 


4900 


M. (z=0) lexp(-z/H, )dz=H. M, (z=0) (1-exp(-4900/H; )) (4.72) 
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the atmosphere is simply given by 


1-R, =1-exp(-4900/H, ) eee, 


Some numbers for the ratio R; are 


3500 4000 


H,(m) [1000 1200 1400 2000 2500 3000 7100 
Rie as Plog UN eRGaz wl G20 bRONOSmaGni4pban7GnrOr25) Geom oWso 


In Chapter 2, H, for aerosol was found to be = approximately 
1200m on a worldwide scale and 1400m in Alaska. A value of 
R, ranging at around 0.03 should thus describe the situation 
sufficiently. 

The value of H,, for water vapor was empirically related 


by Reitan (1963) to the surface dewpoint temperature, t,, as 


following: 


H,,(m) 2300 2400 2600 2900 


For the cases considered at Resolute the surface dewpoint 
temperature is on the average at around -25°C which gives H, 
=3500m by extrapolation and thus R,=0.25. 

It is obvious that the ratios R; for aerosol and water 
vapor are only rough estimates. In the consideration of 
errors at the end of this Chapter this will be taken account 


of. But it can be seen that the previous assumption of 
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constant mixing ratios throughout the atmosphere is rather 


unrealistic. 


4.6.2 Volume extinction coefficients 

The volume extinction coefficients are defined in a way 
similar to the one-layer atmosphere. Coefficients for the 
upper layer (primed) and the lower layer (unprimed) have to 
be distinguished. w is the precipitable water R, and R, mean 
the fraction of aerosol and water vapor, respectively, in 
the upper half of the atmosphere. 

To understand the transmission functions T, and T, it 
must be remembered that the air mass m acts essentially as a 
magnification factor as far as water vapor transmissivity is 
concerned. To calculate water vapor transmissivity for one 
half of the atmosphere’ the precipitable water content of 
this half has to be multiplied by m, and not by m/2. In case 
of Rayleigh scattering, however, which is a function of m 
only, m/2 rather than m has to be taken to calculate the 
transmissivity of one half of the atmosphere. 

The aerosol volume extinction coefficients for the 


upper and the lower layer are: 
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(a= lIeRe te (4.75) 
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The water vapor volume absorption coefficients for the upper 
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The water vapor volume absorption coefficients for the lower 


layer are 
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The Rayleigh volume scattering coefficients for the 


upper layer are 
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The Rayleigh volume scattering coefficients for the 


lower layer are 


Pro-aein (4.91) 
By=4.56*10 (273/T) (p/100) (4.92) 
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The results of the two-layer two-stream approximations 
are given in the next section. The diffusivity factors were 
taken to be 1.66 for reasons which will be explained 


shortly. 
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4.7 Results 

The calculated medians of w, for 288 cases of cloudless 
skies and all other relevant data available at Resolute 
calculated with the Davies-Hay model, the Kondratyev-model, 
the 1.66-model, and the two-layer model are presented in 
Table 8. 

Table 8 shows also the limits within which 68% of all 
solutions are found. Numerical solutions for W, were 
searched in the interval (-1,4). But solutions could not be 
found in all 288 cases. The number of solutions found in 
this interval is also given in Table 8. In some instances, 
solutions were found outside the interval (0,1) within which 


Ww has to be located by definition. The percentage of 


solutions outside this interval, i.e., the percentage of 
physically unacceptable solutions, is listed in the 
following line of Table 8. Finally, the correlation 


coefficient between w, and cos@ is presented. This 
correlation coefficient should be close to zero because any 
deviation from zero is more_ probably due to some faulty 
zenith angle dependence of the transmissivities or the 
related volume extinction coefficients rather than to any 
real zenith angle dependence. 

The Davies-Hay model is obviously insufficient, as it 
yields 42% of all solutions outside the interval which is 
physically meaningful. The scatter of the solutions around 
the median of 0.95 is also very great and the correlation of 


-0.19 indicates that %, is over-estimated for low zenith 


eaeftuoln to sees OBS 107 wo) Fo er 
 ‘weuloseh ts sfaetievs efab feet * | 
[ébon-veyisqanos ed? , haber ra 


Ht pei negeta ate fsbon neyel owt bates 


aid 


an 
a iT a. --) 
foe " 


a Ber) 
- 8 pas 
iia to, x68 etm nF AD tw dine, ert? cals aware - 8 slaer 
anew A 10% e@hdtivioe pear aol vbrwet hel ont 
ad Jon Bikes emartubae tue ha, | Tyrese At he 


di bnyot andttuloe yo. seein anit ~@eeso 88S is mea 
sparnaleni smae! MM -/3 eidet ni navip oats wet fewredrit etds 


i 
’ 


dotrw nietiw (1.0) fevvedrrt caelt abletwa buat? stew eroraute ; 


40 wpednesred SAT neti inttey Ne ‘betecor ed oF ese re 
% sasineoreq et? ..9,4 fisvaagab abt ebt2tue: anotiuloe 
wi nf batetl et ,enptautee stdsigeaseru vi Faotargel 


aortela vis ert vitentiy sical +o enti pai: 
f ; -_ 


atc? .betnezesq et Seon DAR feewied netott om 
yes oausced o18s oF seold od) Biwerde erat teos ies 
 ¥Pfuet amos oF eub yidsderg, @%om at ones mrt nottsts 
aft 10 astiivteaimansad | arid, to sonebneqe) of pias 
yng of nant edie atnetot eco Aor tof lage amy tov bels 

.sonsbreges signs wine 
$} a8 Inetottwent vievotvds et letcm yet-estved oo 
2) fotdw favsaini oft eblietuo enot tutes aeatias | oe 

brusns anotiulae sry) hate oce ——- | 


' 


95 


Davies—Hay two-stream approximation 


one layer two layers 


model 


median and 
68%-—limits 


number of solutions 
in (— 1,4) 


% of solutions 
outside (0, 1) 


Se Po aa ar a 


Table 8: Results of w. at Resolute using four different radiation transfer models. (R,: 
fraction of aerosol in upper half of the atmosphere; R,: fraction of water vapor in upper 
half of the atmosphere) 
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angles. It is , in fact, the region of low zenith angles 
which is responsible for most of the solutions well above 
unity, indicating that the single-scattering approach loses 
its validity for the sun being close to the horizon. This 
result was expected. 

The two one-layer two-stream models (Kondratyev-model 
and 1.66-model) give far more realistic solutions, as only 
as little as about 15% of the solutions are _ physically 
unacceptable, and the scatter of the points is considerably 
less than in the Davies-Hay model. The only major difference 
between these two models is the cos@ -correlation 
coefficient: The Kondratyev-model gives -0.38, which makes 
it unacceptable, whereas the 1.66-model gives 0.06; thus it 
is almost independent of the zenith angle. With respect to 
these criteria the 1.66-model gives the better results. 

The two-layer model for R,=0.03 and R, =0.25 has’ the 
advantage that the scatter of the solutions is lesser than 
in the one-layer two-stream models, but solutions can be 
found only for 64% of all cases from which 23% are 
physically unacceptable. The cos@-correlation coefficient 
has slightly increased to 0.10. 

The fact that in the latter model no solutions can be 
found for one out of three cases needs more consideration: 
The lack of any solution in the interval (-1,4) indicates 
physically that, whatever single-scattering albedo is tried, 
the observed radiation field cannot be matched with the 


theoretically predicted one. piowrersu, tt Inrignt™ be 
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conjectured that complex solutions for the radiation fields 
DT(z=0) and Dd(z=0) should also be considered. So let us 


examine equation 4.18. r becomes complex if 


(a,-a,,) <4(a,a,-a, ay) (4.96) 


For all two-stream models except the Kondratyev-model there 
i L 


is only one diffusivity factor Mmb=nw Sandkethus!B, r=Bai.0S This 


leads to a second-order inequality in w, with the solutions 
a” 

Wy >| - AY =a + VR +Pe )* -4 PRBe+4 (Be-BY)*)/48,(1-Ba) (4.97) 

uc = ( B+ Pe VO pr+ BD? 475th a (AE BE) )/4B.(1-B.) (4.98) 


Examination of expression 4.98 shows that this inequality 
can never be fulfilled as the numerator is always negative 
but w, must always be positive. But expression 4.97 may have 


a positive numerator. Actually it is positive for 


(oh ah)" > Px Pe (4.99) 


In Section 4.4.2 the water vapor volume absorption 


coefficients were introduced as quantities of the order of 


oy BS Psd iy 
magnitude of 10-10 m. As pw and fiw are similar in size 


‘and the square of it 


m_. Ri was introduced in Section 4.4.3 as 4.56410 m, 3k 


their difference is typically 10°m 
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10° 
. <5 go le -¢ -2 
is typically 2-3*10 m , so /3e/3f is around 10m. Therefore, 


ale 


au tel 8 pet oetcnate work Oem iad 
#t xotomoo) a 


(ap k) ie Se" aaa: - 


Led i t= oghlian re: 

evarit eeeuameias ert} iqenie: alebam msente-ow? ft 
BYAT a, ce iF aunt orm dine tn ofan yd bevaut? th 
smotiwlos etl wiftw .@. rt vit taupe directness a 
et eat 


mead tert aan 
(Be.R) (abet WO RES +\+Haae 

ae 
si}teueen! sth tery awode @0) nateaerqpe tojdekiadin 
bvitepec evew!ls sf 1) stem ont es ball itMut ed “ever nso 


evel yver'TO.4 noteestgns FUE -avitieog sa wyewl 6 - 
10? evidteon ef ff 40 \autot. d0te7emun ave 


Lee. e) Py \s (RB 


este nt af inte ote —— , 

$). Yo. ermupe aii + sahenca 
ey eaten - Ei hdog 
aaroreredt im OF tht 


e]s) 


the neglect of complex solutions, although they are 
theoretically possible, does certainly not explain the lack 
OFM3be) Of all So lut tons: 

Another possible explanation is’ the distribution of 
aerosol and water vapor in the atmosphere. Calculations have 
been made for various combinations of R, and R,, the 
fractions of aerosol and water vapor, respectively, in_ the 
upper half of the atmosphere. Table 9 shows the results: The 
upper line in each square indicates the median and _ the 
limits within which 68% of all solutions can be found. The 
second line shows the number of solutions (288 would be the 
maximum possible) and the third line shows the percentage of 
solutions which are outside the- physically meaningful 
intervals c0-s 

Three features are interesting to note: 

Firstly, the median of w, for .R, = Rye Oe oUetss 0,60.) ihe 
one-layer 1.66-model which has inherently the same aerosol 
and water vapor distribution gives w, of 0.80. The 
difference is due to the newly defined volume extinction 
coefficients. Sample calculations with cases of extremely 
low precipitable water show that approximately one third of 
the difference of 0.11 is due to the new Rayleigh 
coefficients the rest must be attributed to the new water 
vapor coefficients. 

Secondly, variation of R, with R, being constant does 
hardly influence the number of solutions. The medians are 


only slightly lower for higher R, and vice versa. 
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“| 0.93(+0.11,-0.16) 
175 
26 


O30(-OM1 - O75 
184 
18 


-0.96(+0.10,-0.17 
178 
33 


0.93(+0.11,-0.16). 
185 
33 


0.89(+0.14,-0.16) 
213 - 
22 


0.88(+0.16,-0.17) 
246 
22 


O76 5(-0:1.75-0416) 
247 


WZ 


0.83(+0.17,-0.16) 
262 


OS 3(+017/-07-7) 
261 
14 


0.79(+0.17,-0.16) 
261 
13 


0.7 56(+0 16°20 J 7) 
271 
11 


0.74(+0.18,—0. 16) 
269 
8 


0.69(+0.20,-0.17) 
375 


0.73(+0.20,-0.18) 
277 
g 


0.64(+0.2 1,-0.19) 
280 
5 


O55(-022,-0139) 
279 
2 


Table 9: Results of “, from the two-layer two-stream approximation with various 
fractions of aerosol (R,) and water vapor (R,) in the upper half of the atmosphere. The first 
line in each square gives the median and the range within which 68% of the solutions can be 
found. The second line gives the number of solutions in the interval (— 1,4) and the third line 
the percentage of solutions outside the interval (O, 1). 
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Thirdly, and probably most important, the number of 
solutions increases for R, being constant and R, increasing 
from 0.03 to higher values. The number of solutions reaches 
its maximum at the rather unrealistic value of R,=0.75 but 
even for R,=0.30 261 solutions can be found. This finding is 
aeesitrongd sindication that “kv=0..03) 1s too low for the Arctic. 
The question that arises now is which value should be taken 
for R, . The more R, exceeds 0.03 the more it deviates from 
the commonly accepted vertical aerosol distribution, but = so 
much the more it approaches the situation for which all 288 
solutions can be found or, in other words, 288 radiation 
fields can be explained theoretically. 

At this point we may remember the problem of 
calculating the aerosol optical depth from integrating 
nephelometer measurements of the aerosol volume scattering 
coefficient (equation 2.14). It was concluded that the 
aerosol optical depth was low by approximately half an order 
of magnitude (which corresponds to a factor of 3). This 
could have been remedied by introducing H,=3600m rather than 
H, =1200m in equation 2.12 which was not done because no 
foundation could be found in literature. 

Assuming now that Shaw’s (1975) vertical aerosol 
distribution for Alaska is not representative of Resolute, 
and assuming H,=3600m instead of 1200m we find from equation 
4.73 a value of R, of 0.26 rather than 0.03. R, =0.26 comes 
closest to R, =0.30 in Table 9 which yields already 261 


solutions (=91% of the total of 288 cases) from which only 
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14% are physically meaningless. Therefore, it may be 
concluded that R,=0.30 comes closer to reality than R,=0.03. 
For the rest of this Chapter we want to accept R, =0.30 and 
thus a value of 0.81 as the median of the single-scattering 
albedo. This result has a larger scatter than the previous 
result of .0.93 (for R, =0.03) and even a correlation 
coefficient of 0.22 (see Table 8). The hourly results’ are 
listed in appendix A. 

It is interesting to see that this new result is very 
similar to the result of the one-layer 1.66-model which 
yielded @=0.80 with a slightly broader scatter but more 
solutions (276 instead of 261) and a considerably smaller 
Cos@ correlation 800 .06.ssinsteadi;s Ofte, 9022) omenthuSied the 
introduction of a two-layer model does not give better 
results than a one-layer model with constant mixing-ratios 
for aerosol and water vapor, as long as the real values 


for R, and R, are not Known. 


4.8 Errors and preliminary discussion 

For the rest of this Chapter it shall be examined to 
what degree errors in the variables change the final result 
Ofeaio=0.9 1% 

If R, deviates from 0.30 by as far as +0.20 the final 
result will change by #40.08, for a deviation of +0.10 it 
will change by 70.05 (Table 9). If R, deviates from its 


assumed value of 0.25 by +0.10 the final result will change 
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by only 0.02. The accurate Knowledge of R, is not as 
important as the precise Knowledge of R,. 

Table 10 shows how the final result changes if the 
following errors are assumed for the variables listed below: 
-radiation measurements (+5%) 

-solar constant (+1%) 

“diffusivity factors »(.1.50v0r 2/00, instead of 1.66) 

-water vapor volume absorption coefficients(+20%) 

-Rayleigh volume scattering coefficients (+10%) 

-aerosol volume extinction coefficients (+10%) 

-downward and upward to total aerosol scatter ratios (+10%) 
-scale height of the atmosphere (+5%) 

Beside the 5%-error in radiation measurements and the 
1%-error in the solar constant which have a foundation in 
literature, as mentioned in Chapter os, similar 
justifications for the errors of the remaining variables 
cannot be given. 

As the Kondratyev-model gave bad results compared to 
the 1.66-model the diffusivity factors for isotropic 
radiation are probably better than those used in_ the 
Kondratyev-model. As was mentioned in Section 4.4.1, these 
isotropic diffusivity factors range between 1.2 and 2.0. But 
inspection of Kondratyev’s diagram (1969,p.19) mentioned 
earlier shows that limits of 1.5 and 2.0 should under’ any 
circumstances be sufficient for the present error limit 


determination. 
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error of | resultant error of | resultant 
variable | change in We| variable | change in o, 


measured values of radiation fluxes 
solar constant 


mi= ml=2.00 
mt= m= 1.50 


all water vapor volume absorption coeff. 
all Rayleigh volume scattering coeff. 
all aerosol volume extinction coeff. 


downward and upward to total aerosol 
scatter ratios 


scale height of the atmosphere 


all the above errors simultaneously 


Table 10: Changes in the final result of w, if errors are introduced separately for each 
variable and then for ail variables simultaneously. 
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Most difficult to justify is the choice of 20% error 
for the water vapor volume absorption coefficients and of 
10% error for the Rayleigh volume scattering coefficients 
and the aeroso| vo lume extinction coefficients, 
respectively. A very simple estimation may show that these 
limits are at least not too low: Let us calculate fy and (3p 
for the case m=2.9, w=2.3mm, R,=0.25, and H=7100m, and _ for 
model atmospheres with increasing numbers of layers. Let us, 
for example, look at the layer which contains or which lies 
immediately above the 75kPa surface. Then we find for /3w and 


Ba in an n-layer atmosphere 
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Application of a two-layer atmosphere’ rather than a 


oO 


multi-layer atmosphere introduces 8% error’ in P. and 1% 
error in /3, for the layer immediately above 75kPa. 
The same calculation for the layer containing or lying 


immediately above 25kPa yields 
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DeOw 


Cooemes tes6" 13.9201 4523 


In this case application of a two-layer atmosphere gives 64% 
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error in (3, and 2% error in fe. 

Excessive errors’ in (e occur only for  downwelling 
radiation in the upper half of the atmosphere, due to. the 
fact that even tiny amounts of water vapor’ induce 
considerable absorption. In a multi-layer atmosphere /3, would 
thus rapidly drop from very high values at the top (0kPa) to 
relatively small values at 50kKPa. Therefore, the average 
value used in_ the two-layer atmosphere underestimates the 
real conditions in a very thin layer at the top, but 
overestimates throughout the rest of the upper half. 

As upwelling radiation has already traversed a_ long 
path before it gets reflected at the ground, considerable 
errors in 3, can be excluded. As only two of the six water 
vapor volume absorption coefficients introduced in equations 
4.76 to 4.85 exhibit large errors, and the remaining four 
are pretty accurate, a 20% overall-error for the Pw- 
coefficients may be sufficient. 

Although the errors in /3, are only about 1-2%, an error 
of 10% was assumed because the [3p -factors calculated for 
diffuse radiation according to equation 4.43 depend on the 
spectral flux data given by Kondratyev which, by no means, 
have to have general validity. 

A 10%-error in fa was assumed because Be is directly 


q 


derived from which had error bars of approximately 10% in 
Chapter 3. 
A 10%-error for the downward and upward to total 


scatter ratios corresponds to the difference between various 
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literature values of this quantity as outlined in Section 
4.2. 

A 5%-error in the scale height H of the atmosphere 
corresponds to an approximate ernor.wof ache inwithe 
determination of an atmospheric mean temperature uf 

Table 10 shows that the source of greatest error is the 
uncertainty in radiation measurements. If the computer 
program is run with all the variables altered in such a way 
that the errors are additive then a total error of 
+0.26/-0.21 is the result. As, by definition, w,<1.00 the 
final result for R,=0.30 and R_=0.25 is 


(JgeQr. Bh HAO tS ta0'. 24) (4.100) 


As R, and Ry, are based on estimations rather than on 
calculations another error has to be investigated. 
Inspection of Table 9 shows that w, increases for decreasing 
R. and decreasing R, but decreases for increasing R, and 
increasing R,. As the value of 0.30 accepted for R, is 
already 10 times greater than the value which has a secure 
foundation in literature it may be sufficient to consider 
only values less than 0.30 for the determination of the 
total error in w,. But as the upper error limit of w, has 
already hit the physically meaningful limit, namely 1.00, 
any reduction of R, from 0.30 to lower values does not 


influence the error limits. 
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The value of 0.25 accepted for R, is again more of a 
maximum possible value than of a mean value for which errors 
have to be considered into both directions. Going back to 
Chapter 4.6.1 we extrapolated the water vapor scale height 
H, far beyond the values given by Reitan (1963) and, as a 
result of that, gained H,=3500m and R, =0.25. Again, lower 
Values? of = Ry do not™ int uence the” upper error Vimi't- as it 
has already reached 1.00. 

We can thus conclude that the aerosol single-scattering 
albedo at Resolute is 0.81. But under the worst conditions, 
especially under the condition of radiation measurements 
being as much wrong as 5%, the single-scattering albedo may 
take any value greater than 0.60. 

With respect to Table 5 an aerosol single-scattering 
albedo of 0.81 is lower than what can be expected for remote 
sites but it is higher than what is found close _ to 
industrial sites. Any interpretation beyond that cannot be 
justified with regard to the large error limits. 

Unfortunately, the cloudless hours at Resolute are 
almost exclusively restricted to spring, so a between-season 


comparison cannot be made. The results given by months are: 


4/78 5/78 3/79 4/79 4/80 5/80 
267 B13 


Orod 0.84 O02 7.0 0.84 0.80 0.74 


The only striking feature in this table is’ the sharp 
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increase from March to April 1979. But this is probably due 
to the cos §-dependence of the result rather than to any real 
phenomenon, because the cos0-mean value for March 1979 is 
Orem buile ToOcwApn?) 1.197 Suit sis. 0i34- 

Further discussion of this result, especially with 
respect to the question of the origin of the arctic aerosol, 
has to await the result of the calculation of the aerosol 


optical depth at the other four stations. 


a 


sd xidedatq et ait nae 


fees vow of narit nwiier, at 
ar O38! Asa 4Ot euti BV 


Attw Wietsegee Shes”. Se Aes 19 


losers otfo ws 847. Yo ofging at soca at 
focoqwn ett toonGttelvohae’ eee v9 ttoess att at aw 


srotceta qua? a caialiaadiaail 


‘6 mie 7 


7) a 


5. Aerosol optical depth at four other arctic stations 


Spiel CRhOGUCTE.TOn 

So far, the aerosol optical depth and the aerosol 
single-scattering albedo have only been calculated for 
Resolute. Resolute was chosen because it is the only 
Canadian station north of 60°N with regular measurements of 
the global solar, the diffuse solar, and the reflected solar 
radiation. 

To get a better understanding for the areal and 
seasona | aerosol turbidity pattern the aerosol 
characteristics should be determined for other’ stations. 
There are thirteen more stations in the Northwest 
Territories where global solar radiation is being measured. 
Data for four of them were obtained for’ further 
calculations, namely for Alert (82.5°N, 62.3°W) which is the 
northernmost meteorological station of North America and 
only some 800km away from the North Pole; for Inuvik 
(68.3°N, 133.5°W) which is the station closest to Alaska and 
thus most readily comparable to Barrow, Alaska, where most 
of the Arctic aerosol research has been centered to date; 
for Baker Lake (64.3°N, 96.0°W) one of the southernmost 
stations of the Northwest Territories and one of the few 
continental ones, and finally for Frobisher Bay (63.7°N, 
68.6°W) which is located in the southeasternmost portion of 


the Northwest Territories and thus closer to the Atlantic 
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Ocean than any other arctic station (see Figure 1). 

Because only one, rather than three, radiation fluxes 
are measured at these stations two more variables have to be 
inserted as constants in the governing equations. The lack 
of reflected solar radiation measurements suggests that the 
albedo has to be estimated. The lack of diffuse solar 
radiation measurements necessitates that one of the two 
variables, aerosol optical depth or aerosol single- 
scattering albedo, has to be taken as a Known constant for 
the calculations. 

From what has been said about these variables it is 
obvious that the aerosol single-scattering albedo must be 
taken as a constant, so the models can be solved for the 
aerosol optical depth. Even at this point a warning must be 
given that, due to the introduction of two more constants, 
the results of T, for the four other arctic stations are 
inevitably less accurate than the results of . for 


Resolute. 


5.2 Radiation models 

Two radiation models are examined for the purpose of 
determining the aerosol optical depth €,: the Davies-Hay 
mode | combining the equations for direct and diffuse 
radiation and the one-layer two-stream approximation with 
diffusivity factors of 1.66 and a direct radiation term. As 


both models describe the radiation fluxes for cloudless 
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skies, only hours with 10/10 sunshine duration and 0/10 
cloudiness can be used. 

For the aerosol single-scattering albedo w~, a constant 
value of 0.95 was used in the Davies-Hay model and a 
constant value of 0.80 in the 1.66-model. These values 
correspond to the results of w, gained for Resolute with the 
Davies-Hay model and the 1.66-model, respectively. 

The appropriate choice of surface albedos needs more 
extensive discussion: The surface albedo is to a certain 
degree dependent on _ the cloudiness of the sky (see, e.g., 
Wallén, 1948; Maykut amd Church, 1973; Petzold, 1977). 
Therefore, only cloudless hours were taken to determine the 
albedo at Resolute. As all cloudless hours at Resolute 
coincide with complete snow cover, only the albedo for the 
snow cover period could be calculated and was found to be 
Oar2s Ata-this: rpointori termus tiuebe remembered that the 
radiometer is mounted close above the ground. Therefore, 
this result is only representative of the immediate vicinity 
of the radiometer. Thus, choosing albedo values 
representative of the entire Arctic for all seasons requires 
a brief review of the vast albedo literature. 

Beside the usual textbook values for the albedo of a 
tundra surface in summer and in winter, the following 
investigations were found for the North American Arctic: 

Jackson (1961) reports on helicopter measurements made 
in Labrador-Ungava. For areas covered by bare rock he found 


an albedo of 0.08 to 0.12 in summer and 0.70 in winter. 
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Davies (1965) reports that in the same region a_ value 
of 0.15 could be found for the albedo of tundra surfaces in 
summer. He made also aircraft measurements. 

Ground-based measurements at Moosonee (northern 
Ontario) gave values of 0.17 for summer and 0.72 for winter 
(Moller, 1965). 

Aircraft measurements were made by McF adden and 
Ragotzkie (1967) over northern Canada. The following results 
are valid for clear skies: snowcovered tundra 0.89; tundra 
without snow but with partly frozen lakes 0.25, ...unfrozen 
lakes 0.11, ...nmo lakes at all 0.15. 

Ground-based measurements in the vicinity of Baker Lake 
made by Ahrnsbrak (1968) give summer values of the albedo of 
0.14 to 0.24. 

Wilson and MacFarlane (1969) found that the albedo at 
Porte-de-la-Baleine (northern Quebec) dropped from 0.85 
after a snowfall to 0.77 two weeks later. Measurements were 
made on ground. 

Finally, Maykut and Church (1973) in a study made about 
Barrow, Alaska, give winter albedos of 0.80 to 0.90 
(average: 0.84) and a summer average of 0.18. Measurements 
were also made on ground. 

With these results in mind and Knowing that all four 
stations are in a tundra environment (Canada, Dept. of 
Energy etc., 1974, p.45f.) it was decided to use 0.80 as the 
albedo for hours with snowcover, 0.15 for the period without 


snow, and to omit all occurrences with traces of snow or 
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unknown snow conditions. 

Let us now return to the radiation models: One more 
difficulty arises from the fact that ozone soundings are not 
available for these stations. Monthly mean values’ for 
Resolute for the years 1975, 76, 78, 79, and 1980 were 
calculated from published ozone data (Canada, Environment 


Canada etc., 1975ff.) and graphically smoothed to give: 


430 465 470 465 430 385 350 320 300 310 345 390 


These mean values compared with daily ozone readings 
from the ozone stations at Barrow (Alaska), Churchill 
(Manitoba), and Goose Bay (Newfoundland) indicate that in 
almost all cases deviations of the real ozone amounts’ from 
the amounts given above are less than 20%. 

Now all necessary data are available for use in_ the 
above-mentioned radiation models. 

Addition of equations 3.3, 4.2, 4.3, and 4.4 yields the 
global solar radiation flux on the ground according to the 
Davies-Hay model. This equation can immediately be solved 
for T, and thus for t,. As the equation is of first order 
there is always one solution except for the cases of 
negative T, which are physically meaningless. 

In the one layer two-stream approximation the global 
solar radiation flux on the ground is given by equations 3.3 


and 4.21. To solve the sum of these equations numerically 
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for 8, the relation T, =H, has to be substituted into 
equation 3.5 and the latter into equation 3.3. Furthermore, 
as the splitting of the global flux into the direct and 
diffuse fluxes is not Known, a ratio of 2:1 was assumed 
(which turns out to be quite accurate for many cases at 
Resolute) to calculate the Rayleigh volume scattering 
coefficient for upwelling diffuse radiation according to 
equation 4.49 and the water vapor volume absorption 
coefficient for upwelling diffuse radiation according to 
equation 4.53. Numerical solutions FOme og were searched in 
the interval (-1*10%, 2*10'm.) which, multiplied by the 


atmospheric scale height H, give the aerosol optical depths 


5.3 Results 
Both models were applied to all four stations and the 
following statistical parameters were calculated (Table 11): 
-the correlation coefficient between 1, and the cosine of 
the solar zenith angle, cos@, 
-the intercept, 3%, and the slope, b, of the regression 


equation 
T,=a+b*cos@ (Sali) 


-the median of ©, and the range within which 68% of all 


results of ©, can be found. 
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station model! correlation median and 68% limits of 
T,# COs 0 Ce 


0.100 (+0.054,-0.056) 
0.092 (+0.034,-0.034) 
0.087 (+0.074,-0.042) 
Ostety7i(40:0507-0.052) 
0.078 (+0.056,-0.058) 
0.082 (+0.06 1,-0.042) 
0.098 (+0.056,-0.045) 
0.102 (+0.045,-0.042) 


Davies—Hay model Ol125 
0.089 
O23 0.068 | 0.093 
1.66-model 0.62 0.05242 08162 
Davies—Hay model 0.00 0.082 | -0.001 
Baker L. (324) 
1.66-model 0.44 0.047 | 0.121 


Davies—Hay model 


Alert (303) 


1.66-model 


Davies—Hay model 


Inuvik (745) 


Frobi. B. (394) 


1.66-—model 


Table 11: The correlation coefficient as well as the regression constants, a and b, of the 
regression T. vs. cos@, and the medians and 68% limits for Ta calculated with two 
different radiation models at four arctic stations. The number in brackets behind each 
station name indicates the number of observations. 
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Let us first look at the correlation between 7, and 
cos§@. The Davies-Hay model yields small correlations whereas 
the 1.66-model gives correlations of up to 0.62. It is 
conceivable that this high positive correlation represents a 
seasonal march of ‘7%, with high values in summer and low 
values in winter. But inspection of the data shows that a 
strong cos@ -dependence of T, can also be found on a daily 


co 


basis. As it is very unlikely that %,. increases in the 
forenoon and decreases in the afternoon the cos0-dependence 
or at least a big portion of it is due to a weakness of the 
radiation model. So the simpler Davies-Hay model gives more 
reliable results. 

If it is assumed that even the cos@-dependence of the 
Davies-Hay model is faulty then the maximum error due_ to 
this dependence can be estimated using the regression 
coefficients: The lowest monthly mean values of cosO are 
found to be 0.1 and the highest ones, e.g. in Inuvik, are 
around 0.5. For Inuvik the monthly T, mean values would then 
be’ 0.068+0.1*0.093 = 0.077 and 0.068+0.5*0.093 = 0.115, 
respectively. In other words, under the worst circumstances, 
an amount of 0.115-0.077 = 0.038 of the difference of T, 
between summer and winter months would be due to. the 
erroneous cos§-dependence of the final result. But, as will 
be seen shortly, the differences betwen months are in-= many 
cases much greater than 0.038 and in some cases the low Jl, 
mean values can be found together with high cos@ mean values 


and vice versa. Furthermore, inspection of Table 11 shows 
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that Baker Lake has no cos §-dependence at all, and Alert has 
even a negative cos@-dependence. If the Davies-Hay model 
really contained an erroneous cos@-dependence it should be 
found at all stations. Therefore, it may be concluded _ that 
the Davies-Hay model is largely free of any systematic error 
due to the solar zenith position. 

Table 11 shows also that the medians gained for two 
stations using the two different models are very similar, 
namely for Baker Lake and Frobisher Bay. The two numbers in 
brackets give the range within which 68% of all results can 
be found. The Davies-Hay model gives slightly larger ranges. 
This need not necessarily be a disadvantage as _ certain 
seasonal changes in CT, should be expected. 

For the rest of this Chapter the results calculated 
with the Davies-Hay model will be used. 

The Davies-Hay model with fixed values for the surface 
albedo and the aerosol single-scattering albedo was used to 
recalculate T. for Resolute. This procedure yielded a 
correlation between T, and cos@ of 0.27 as compared to 0.24 
from the direct radiation calculations of Chapter 3. The 
median and the 68%-limits for all cloudless hours amount to 
0.103 (+0.095, -0.062) as compared to 0.087 (+0.052, -0.030) 
from the direct radiation model. It seems that’ the 
Davies-Hay model slightly overestimates T,. This should be 
kept in mind to compare the results from the four global 


radiation stations with the result from Resolute. 
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Tables 12 to 15 show the monthly medians for T, 
calculated with the Davies-Hay model! for Alert, Inuvik, 
Baker Lake, and Frobisher Bay, respectively. The line under 
each median indicates the 68%-range. Results are given only 
for months with at least ten single results or with results 
stemming from at least three different days. The tables wil] 
be commented in Section 5.5. Appendix B contains the hourly 


input and output data for these stations. 


5.4 Errors in the results 

The Davies-Hay model contains several variables with 
considerable inaccuracies. They have to be subjected to an 
‘error analysis. 

First of all, an error of 5% has to be assumed for the 
global radiation measurements and an error of 1% for the 
‘solar constant (see Chapters 3 and 4). 

An error of 10% is assumed for the upward and downward 
to total aerosol scatter ratios as in Chapter 4. 

For the ozone content an error of 20% was_ taken 
corresponding to the findings outlined in Section 5.2. 

For the surface albedo the lower error limit was’ taken 
to be 0.05 in summer and 0.65 in winter, and the upper error 
limit was taken as being 0.30 in summer and 0.90 in winter. 
All albedo investigations mentioned in Section 5.2 fall into 


these limits. 
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Table 12: Monthly medians of T, at Alert The numbers in brackets indicate the 68% range 
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Table 13: Monthly medians of T, at Inuvik. The numbers in brackets indicate the 68% range 
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Table 14: Monthly medians of T, at Baker Lake. The numbers in brackets indicate the 68% 
range of T,. 
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Table 15: Monthly medians of ©, at Frobisher Bay. The numbers in brackets indicate the 
68% range of T,. 
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Finally, for the single-scattering albedo the limits 
were taken as 0.80 and 1.00. It must be remembered that the 
Davies-Hay model gave a value of 0.95 for w, at Resolute, so 
0.80 represents a considerable negative deviation. A value 
of 1.00 is the upper limit which is meaningful. 

A justification that the errors of the other variables 
in the Davies-Hay model are negligible was given earlier. 

Table 16 shows how the median of 7, for Alert, which is 
0.100, changes if these errors are introduced separately. 
The lower part of Table 16 shows what happens to the medians 
of all four stations after all errors are introduced 
simultaneously so that they are additive. It turns out that 
these total errors are huge. The negative error deviations 
have approximately the sizes of the medians themselves (see 
Table 11) and the positive deviations have even about twice 
the sizes of the medians. In other words, if the median of 
T for a station is called m then the total error analysis 
suggests that the real median is located somewhere between 0 
and 3m. As the medians are around 0.1 this result means that 
the error bars’ have to be extended from about 7. = 0.0 to 
SGD te ee at) aoe Wh Ch makes the results completely 
meaningless. 

But inspection of Table 11 shows that the medians of 
the four stations are very similar, namely 0.100, 0.087, 
0.078, and 0.098, respectively. If the error bars were as 
large as indicated in the last paragraph it would be very 


unlikely to get results of such close similarity. Table 16 
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error of resultant error of resultant 
variable | change in7,| variable | change int. 
measured values of global radiation 
solar constant 


downward and upward to total aerosol 
scatter ratios 


ozone content 
surface albedo in summer 0.05/0.30 
surface albedo in winter 0.65/0.90 


single—scattering albedo 0.80/1.00 


all these errors simultaneously at 


Baker Lake 
Frobisher Bay 


Table 16: Changes in the final result of Ta. if errors are introduced separately for each 
variable at Alert (upper part) and if all errors are introduced simultaneously so that they are 
additive at Alert, Inuvik, Baker Lake, and Frobisher Bay, respectively (lower part). 
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shows that by far the biggest portion of the total error 
originates from the assumption that the error that has to be 
applied for global radiation measurements is as big as _ 5%. 
This error is also responsible for a large fraction of the 
total error of the aerosol single-scattering albedo in 
Chapter 4. 

With a considerable degree of confidence the _ following 
two conclusions may thus be drawn: 

Firstly, a 54 error in the radiation measurements is 
too pessimistic. 


im 


Secondly, the errors in %, are considerable. An 
intercomparison of the absolute @,-values between the 
stations is not justified. But as the errors due to a_ badly 
calibrated pyranometer or due to an over- or underestimation 
of the other variables are systematic, i.e., do usually not 
change sign and extent over time, an intercomparison between 
the seasons at a given station is’ feasible is the 


interseasonal change is considerable and occurs regularly 


each year. 


5.5 Discussion of the results 

The medians for Alert (Table 12) show one pronounced 
feature, namely a relatively high aerosol optical depth in 
the first half of the year and then, approximately in July, 
a sharp drop to values about half as great as before July. 


The maximal aerosol turbidity occurs betwen March and May 
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and its extent is variable from year to year. The low values 
for March 1979 may indicate that the period of high aerosol 
turbidity does not extend all through winter but is confined 
to several months in spring. An influence of the eruption of 
Mt. St. Helens cannot be found for July 1980 suggesting that 
the ‘aerosol cloud’ has not arrived there before the end of 
July. 

The results for Inuvik (Table 13) show a_ similar 
periodic change of aerosol turbidity. As more values are 
available in Inuvik for the months of February and March it 
becomes obvious that the occurrence of high aerosol 
turbidity is confined to a short period of about two months 
in spring. Its onset is variable from year to year but high 
aerosol turbidity may not even occur at all as in 1980. The 
sharp drop from high to low aerosol turbidity occurs earlier 
in the year than in Alert, namely already in May and_ June. 
After an apparent midsummer mimimum aerosol turbidity rises 
again to higher values in late summer and early fall of 1978 
and 1979. From June to September 1980 aerosol turbidity is 
higher than in the previous years which may be due to Mt. 
St. Helens. 

For Baker Lake (Table 14) values were unfortunately 
available for only 14 months. But these data are enough to 
show the same feature that was observed in Inuvik: After 
very low winter aerosol turbidities there is a two month 
maximum in March and April. The second maximum in_ late 


summer is, in the case of Baker Lake, more pronounced than 
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in Inuvik and is at least of the same magnitude as the 
spring maximum. 

Finally, the seasonal pattern at Frobisher Bay (Table 
15) is to a lesser degree periodic than at the other’ three 
locations. 1978 shows a very short but intensive maximum in 
March, a minimum in April and high values’ thereafter. In 
1979 only the short and_ pronounced March maximum can be 
found again. The entire year 1980 shows high aerosol 
turbidities. The May maximum occured before the eruption of 
Mt. St. Helens on May 18th, 1980. 

At this point it should be remembered that the typical 
spring maximum was also found at Resolute (Chapter 3) which 
is thus a common feature for the entire Canadian Arctic. 

In the last Chapter the conclusion was drawn that, due 
to the large errors, an intercomparison of the absolute 
magnitudes of aerosol turbidity between the stations is not 
feasible. But as annual medians were _ found to be very 
similar at all five stations, namely T, = 0.09 + 0.01, this 
value may be compared with the results from other authors 
given in Tables 1 to 3: 

The results of this thesis are in good agreement with 


oo 


other T,-measurements in the Arctic but not with T,-values 
calculated from integrating nephelometer readings in the 
Arctic. The aerosol optical depth in the Canadian Arctic is 
lower than in most parts of the midlatitudes but is 


distinctly higher than in Antarctica. 
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5.6 Interpretation of the results and conclusion 

The purpose of this thesis was to calculate the aerosol 
optical depth and the aerosol single-scattering albedo from 
radiative transfer models at several stations in_ the 
Canadian Arctic in order to contribute to the Knowledge of 
the areal and seasonal extent of the aerosol turbidity and 
to the question of its origin. 

This final discussion will deal with three questions: 
-Are the radiative transfer models used sufficient for this 

purpose? 
-Are there common characteristics of the aerosol turbidity 
in the Canadian Arctic? 
-Where does the arctic aerosol come from? 

As for the radiative transfer models the following must 
be stated: Integrating nephelometer measurements made in the 
surface boundary layer are insufficient for determining the 
seasona | characteristics of the entire atmosphere. 
Therefore, models for the transfer of radiation through the 
entire atmosphere have to be applied. 

The direct solar radiation model used to determine T,. 
at Resolute is indispensible, as no serious. alternative 
exists. The single-scattering model (Davies-Hay model) and 
the multiple-scattering model (two-stream approximation in a 
one-layer atmosphere, or 1.66-model) are not interchangeable 
as the first one has to be taken to calculate tT. from global 
radiation fluxes and only the latter one can be used to 


determine the aerosol single-scattering albedo. In _ other 
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words, Davies and_- Hay (1980) have developed the 
single-scattering model to such a perfection that it is now 
a competitor with the simplest versions of 
multiple-scattering approximations. Application of models 
with more than one layer is not justified as long as_ the 
vertical profile of absorbers and scatterers has to be taken 
from literature values which claim general global 
applicability. 

The major stumbling-block for the application of any 
radiation model, as sophisticated as it may be, is the fact 
that a large error has to be taken into account for’ the 
radiation measurements. But the results of Chapter 5 point 
out that the generally accepted error of 5% is most probably 
too pessimistic. All error bars in this thesis are thus too 
large. It is a necessity to reassess the inaccuracy of the 
radiation measurements per formed by the Atmospheric 
Environment Service with the aim of giving more realistic 
error limits for empirical radiation studies of this kind. 

Even with errors of 5% of the measured radiation’ the 
results can be used to find certain seasonal characteristics 
of the aerosol optical depth. The most outstanding feature 
is a period of high aerosol turbidity in spring which does 
not last longer than three months. In the northern and 
western Arctic the aerosol turbidities in summer are about 
half as great as during the spring maximum. The stations in 
the southern and especially in the southeastern parts of the 


Northwest Territories show a secondary period of high 
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aerosol turbidity in late summer. Due to heavy cloudiness in 
fall and due to the polar night in winter not many results 
could be gained for these seasons. But, from the few data 
available, winter seems to be ae period of low aerosol 
turbidity. This result does not support the assertions of 
various researchers who, by making only few measurements, 
claim that winter is the period of high and _ summer _ the 
period of low aerosol turbidity. 

The results gained in this thesis rather suggest that 
three regimes are interfering: 

Firstly, the background aerosol regime which amounts to 
about €, = 0.05 and is comparable in size to the aerosol 
turbidity of the remotest place on Earth: Antarctica. 

Secondly, there is a spring regime which doubles or 
triples or even quadruples the aerosol turbidity for several 
weeks in spring. The data of 1978 and 1979 show that it 
occurs in the southern and eastern portions of the Northwest 
Territories earlier than it does in the northern and western 
portions. As the Arctic is almost unpopulated and there is 
continuous snowcover throughout the Canadian Arctic’ in 
spring the aerosol must be imported. But it must not be 
concluded that the ‘aerosol cloud’ moves in a southeasterly 
flow over the Arctic as the mean flow would be capable of 
distributing the ‘aerosol cloud’ in much less time than one 
or two months. 

Thirdly, there is a summer regime which affects only 


the southern and eastern parts of the Northwest Territories. 
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This regime is essentially the same regime that is 
responsible for the midlatitudinal summer maximum (see the 
first lines of Table 1). It is due to the fact that any 
sur face free of snow is a potential aerosol source, 
especially in combination with the occurrence of convective 
storms. In the case of the Arctic it may also be due to 
forest fires in the adjacent boreal forests. 

With respect to what was said in Chapter 2 about the 
origin of the arctic aerosol, a plausible explanation for 
the short but intensive spring maximum seems to be the 
duststorm hypothesis. As duststorms do not only occur in the 
semi-deserts and deserts of central Asia it should be 
wor thwhi Le investigating whether dust from similar 
occurrences on the North American prairies can reach the 
Canadian Arctic during spring. 

The . lengthy determination of the single-scattering 
albedo did unfortunately not contribute conclusively to’ the 
solution of the problem of the aerosol’s origin. This is so 
because the final result of 0.81 does not clearly indicate 
whether the aerosol is typically a natural one or typically 
a man-made one. The large error bar of the final result is 
further aggravating its interpretability. 

Furthermore, a warning should be issued at this’ point 
against the use of the aerosol single-scattering albedo as a 
clearcut indicator of the aerosol’s origin. This has’ three 


reasons. 


oy net ; 


et heed. -amtparr ai 
ed ese) mumntaam wean tortie 
vig edt toast -ent os su a 
af;zuou loess lerineroq = Pa 
S20 afd aster ba: 


by aed? Bs Jaw. ++ atrpad sri? io. ait eres 


ri 


iS apne tw 


tear0? | eened inedetbhe ‘ert, rf 

bcc Sash ct bise eaw  seiw od tsauae As W 
‘ave efdlabsiq #., toaGhee offers eri? tov nh yt 

antige:  svhanetnt... oud. dnarlaatel 
1 fon ob emnretarave a eteen ogvn mote 
B1InsS » @lreasth hrs ta 

ru went jae tees t 

- | treat disoh ery fe sino 
tay 2 gniawb. of fowl neth 


sell eth iia at. © porisnimeies yritanst, a na 
Un : a 
Vio Teh ¥ 0d aR NEES | bib, ab: zs 


ae ZIT sniotae @ Woeetes ent, To mph alone and, ‘a not ty 
soibnt \isseld for eeoe Paro to-2fuegn [ett t wine 
ns vilpotay? ef. losers wth ett rs 


— | . ' 
2 r q fiean ok a =}? 77 “ec TTS S78 | anit - gba Ae 
wit! tohetetaneint2ft getiave ws S. 


intoq efdi ts beyea! ot Sfaone got naa & ,ecomnert3 Wy 
. ae 


oy 
* 
a 


6 #5 obedis oni st isos-el pata: Jozo%en, a3 Fe, 
| sawit seri ell .nigiw; 2 ‘foanres ar? to 
fe - | oe 0 
al ae 

oa a a 


~ 
a 


132 


Table 5 showed already that desert dust storms may 
produce relatively low single-scattering albedos due to’ the 
fact that dried and ground remnants of plants’ behave 
optically like man-made aerosols. The same is true for ashes 
from forest fires. Although some efforts have been made to 
determine the imaginary part of the refractive index for 
various substances (see Chapter 2) the same attention should 
also be given to the determination of the single-scattering 
albedo of these substances because the latter can be 
calculated from the former only if the particle 
size-distribution is Known (see Figure 2). It should be kept 
in mind that it is the single-scattering albedo which is 
ultimately the variable used in many radiation and even 
global climatic models. 

The second reason why the aerosol single-scattering 
albedo is not an optimal tracer for aerosols is’ the fact 
that the particle size-distribution undergoes changes with 
the ageing of an aerosol. Mie theory shows that’ the 
scattering and absorption cross-sections of a dielectric 
sphere change rapidly when the radius of _ the sphere 
increases. Therefore, with the change of a particle 
size-distribution a change of the ratio of the _ total 
absorption cross-section to the total scattering 
cross-section per unit mass of aerosol occurs and thus’ the 
single-scattering albedo changes. 

Another reason why the aerosol single-scattering albedo 


might be changeable is its dependence on relative humidity 
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(Hanel, 1976). As this dependence 


relative humidities greater than 80% 


portion of the entire data set used. 


to be considered, too. 
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Appendix A: Data for Resolute 


The columns have the following meaning: 
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line number 

year 

month 

day 

hour 

global solar radiation (kKU/m*hour ) 

diffuse solar radiation (kU/m*hour ) 

reflected solar radiation (KU/m*hour) (****: missing) 


ozone (cm NTP) 


0 oO fA N DO OO FP WD ND 
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: precipitable water (cm liquid) 
11: station pressure (kPa) 

12: station temperature (°C) 

13: cloudiness (tenths) 

14: relative humidity (%) 

15: visibility (km) 

16: wind direction (degrees) (0: calms) 
17: cosine of the solar zenith angle 
18: ozone transmissivity 

19: Rayleigh transmissivity 

20: water vapor absorptivity 

21: aerosol optical depth 


22: aerosol single-scattering albedo 
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Appendix B: Data for the other four stations 
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